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FOREWORD 


This document is Volume I of the final report of NASA Contract 
NASS-24747, "Development of Thermal Stratification and 
De stratification Scaling Concepts, " Mr. L. A. Holmes served as 
the program manager for the initial phase of this study. In 
August 1970 Mr. T. M, Lovrich assumed the responsibilities of 
program manager and is responsible for the performance of this 
study. The study was conducted for the George C. Marshall Space 
Flight Center of the National Aeronautics and Space Administration 
under the technical direction of Mr. T. W. Winstead, Thermo- 
dynamics and Fluid Mechanics Branch, Propulsion Division, 
Propulsion and Vehicle Engineering Laboratory. 

Appreciation is expressed to Dr. S. H. Schwartz, Associate 
Professor Mechanical Engineering, West Virginia Unive: sity for 
his assistance as a consultant. The authors wish to recognize 
Mr. L. A. Holmes, Engineering Instructor, Modesto Jx*. College 
in California, for initiating most of the initial scaling concepts 
analytical work and test design, and to Mr. T. W. Winstead, 
MSFC/NASA Contract Technical Supervisor, for his assistance dur- 
ing the performance of this research-study; and to J. E. Rogan, Heat 
and Mass Transfer Branch Chief, of MDAC for his support in the 
performance of this research study. 

The final report consists of two volumes; 

Volume I Definition of Thermal Stratification Scaling 
Parameters and Experimental Investigation 

Volume II Stratification — Experimental Data 

Volume I is a report of the final results, conclusions, and recom- 
mendations of the study. It includes a presentation of the analytical 
methods used in defining the thermal stratification dimensionless 

Hi 



scaling parameters and the experimental investigation performed 

with saturated Freon 113 in three closed tanks thermodynamically I 

and geometrically scaled for detei'mining the validity of these seal- j 

ing parameters. Volume II contains a presentation of the expei’i- 
mental data collected and utiliaed in this stratification research study. 

n. 

Work performed with the 6 inch and 18 inch diameter tanks was sup- 
ported under this contract, while the work performed with tlie IE inch 
diameter tank was primarily supported by MDAC-West IRAD "Liow-G | 

Fluid Mechanics and Keat Transfer" funding. ! 
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ABSTRACT 


A detailed study to determine the important dimensionless param« 
eters associated with the thermal stratification and pressure history 
of a heated container of liquid and its vapor was made. A test 
program was then designed to verify the validity of these scaling 
parameters. The Modified Grashof number, the Fourier number, 
and an Interface number were the three pi-groups that were 
parameterized using a single test liquid. Freon 113. Three cylin- 
drical test tanks with spherical dome end caps, 6, IZ, and 18 inches 
in diameter with an Li/D of 3, were built. Blanket heaters covered 
the tanks and 77 thermocouples monitored rhe temperatures of the 
liquid, the ullage, the tank walls, and the foam insulation encapsula- 
ting the tank. A centrifuge was used for the 6 inch tank to preserve 
the same scaling parameter values between it and the larger tanks. 
Tests were conducted over a range of Gr'** values, between 10 to 
10^^, and the degree of scaling was checked by comparing the 
dimensionless pressures and temperatures for each scaled pair of 
tests. Results indicate that the bulk liquid temperature, the surface 
temperature of the liquid, and the tank pressure can be scaled with 
the three dimensionless parameters. Some deviation was, however, 
found in the detailed temperature profiles between the scaled pairs 
of te s ts . 
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Section I 
INTRODUCTION 

Many future space missions will rely on long-term cryogenic-chemical pro- 
pellant storage and fuel-transfer systems. The feasibility of these mission 
concepts and design configurations will continue to be strongly influenced by 
the ability to effectively design a cryogenic propellant storage system. Con- 
sequently, it is important that the designer have a strong technical platform 
from which to predict tank design pressures and thermodynamic and fluid 
mechanical behavior in a low-g, near-zero gravity environment. 

Although considerable ground-test data have been gathered, and numerous 
analyses carried out in the area of propellsint thermal stratification (and 
de stratification), the basic scaling parameters that provide the relationship 
between the ground model tests and the full-scale flight vehicle have not been 
adequately substantiated. The associated phenomena are sufficiently com- 
plicated so that useful scaling parameters cannot be established by analytical 
means alone. Thus, a systematic analytical and experimental study of scaling 
parameters for thermal stratification (and de stratification) was made in this 
program. It was the purpose of this program to substantiate the validity of 
the defined scaling groups. 

The conservation equations describing the liquid and its vapor phase were 
nondimen si onalized and the resulting dimensionless parameters were used to 
define a test matrix. The test program was designed to determine whether 
or not the dimensionless pressures and temperatures met the scaling criteria 
for two tanks which have dimensionless parameters equal to one another. 

Three geometrically and thermodynamically scaled test tanks were designed 
and fabricated. A destratification jet orifice was at the bottom of each tank. 
These tanks were capable of heatirig the liquid and its vapor, either separately 
or together at various liquid levels. Seventy-seven thermocouples were 



installed in the liquid, the vapor, the tank walls, and the insulation. Tempera- 
ture and pressure measurements were monitored to verify the scaling 
parameters. 

j 

Difficulties were encountered in obtaining scalable stratification/destratification 
data. These were in the areas of tank material incompatibilities with Freon 113 
and the high-g effects on the instrumentation equipment during testing on the 
centrifuge. Consequently, a new test matrix was defined and the data 
retrieving system was revised midway into the program. The number of 
basic scaling \ests was reduced and a single fluid, Freon 113, was used. 

The new test plan consisted of 20 scaling tests and an additional 17 tests for 
determining data reproducibility. The program tasks were to; (1) define the 
scaling parameters analytically, (2) fabricate the test apparatus, (3) define 
a test program, (4) conduct tests, and (5) reduce and analyze temperature 
and pressure data based upon the defined scaling groups. 

During the tests, the three scaling pi groups, Fo, Gr=5«, and I were maintained 

constant according to the test matrix. These tests covered a range of 

14 16 

Modified Grashof numbers, from 10 ■ to 10 , Fourier numbers, from 

10 ^ to 10 and Interface numbers, from 7. 54 x 10 ^ to 2.42 x 10 
De stratification tests were conducted at the termination of each stratification 
test by activating a pump which produced an axial mixing jet. The jet 
Reynolds numbers ranged from 20, 000 to 45, 100. 

This document discusses the analytical procedures, the experimental equip- 
ment, the test matrix and experimental program definition, testrng, the 
results, and the conclusions and recommendations. 
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Section Z 

DIMENSIONLESS PARAMETERS 


2. 1 INTRODUCTION 

This phase o£ the study was directed towards identification of the dimensionless 
parameters associated with the thermal stratification process. The technique 
of nondimensionaliaing the set of conservation equations and boundary conditions 
which describe the heated liquid- vapor system was used to determine the param- 
eters analytically. This method was selected because it lent itself readily to 
the stratification process and is a rigorous and accurate technique. However, 
it should be noted that the set of dimensionless parameters obtained can only be 
as accurate as the equations which describe the physical process. Further- 
more, it has been assumed that a unique solution to this set of equations exists. 

The first step in the nondimensionaliaation process was to write the conserva- 
tion equations and the boundary conditions describing a specific tank geometry 
and liquid- vapor system. For simplification, the thermodynamic system was 
limited to a flat liquid -vapor interface and only a 2 -phase system was 
considered. The conservation equations were written separately for the liquid 
and the vapor. The boundary conditions include both the tank walls and the 
liquid-vapor interface separating the two phases. 

The next step in the nondimensionalization process was to select a set of 
quantities, characteristics of the physical system, which include all of the 
basic physical units of the system. These characteristic values, such as 
length and time, were then used to nondimensionalize the independent and 
dependent variables included in the governing equations and the boundary con- 
ditions, Selecting the proper characteristic quantities is necessary to provide 
a consistent and useful set of scaling parameters. These values are measurable 
constants that are characteristic, yet distinguishable, and can be related to a 
specific tank. 
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The resuitiug dimensionless variahles are then substituted back into the 

differential equa.tions thus transforming them into a set of dimensionless j 

differential equations and boundary conditions. The dimensionless coefficients 
associated with these equations are the dimensionless parameters. The follow- ; 

ing discussion will cover in detail the steps required to obtain these parameters. | 

- I 

I 

2. 2 CONSERVATION EQUATIONS AND BOUNDARY CONDITIONS | 

A two-dimensional system in cartesian coordinates was used when developing i 

the set of conservation equations and boundary conditions for the thermal | 

'i 

stratification process. However the same group of dimensionless parameters | 

'would be achieved without this simplification. Figure 2-1 illustrates the 

i 

overall two-dimensional system and the boundaries. The governing equations | 

were written in differential form for an arbitrary location both in the liquid 

and the vapor where the gravity vector is in the y direction. It was assumed 

that the transport properties are constant and that liquid is incompressible, ; 


The following set of equations apply to a single component, two-phase system 
with a, flat liquid-vapor interface. Holmes (Ref 1) presented a more general 


FIGURE 2-1 

GENERALIZED STORAGE SYSTEM 
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Set of equations to allow for a curved interface, two components, and 3 
phases. The momentum equations listed below incorporate the Boussinesq 
approximation which gives rise to the gPAT term. In addition, components of 
the total pressure are listed separately. The pressure p at a point is expres- 
sed as 

P = Pref + Ph+Pm 

where Pj.g£ is the pressure at the location where p^, the hydrostatic pressure, 

is zero. The value of p^ is that pressure due to motion and is the dynamic 

pressure. Thus, in the liquid p ^ = p whereas in the vapor p p - 

p„ (i. e. , this is effectively the measured ullage pressure). The gPAT 

^y - Y 

term in the momentum equations comes from the 9pj^/9x term and the body 
force term and only the 9p^/9x term remains. The pressure of the ullage, 
p^, is found in the vapor-energy equation and is defined by the equation of 
state. The measured value of p^ will be close to the total pressure, p, since 
p is small under the conditions of this study. Thus, the values of p^ and 


p^ are entirely sufficient. 


LIQUID REGION : 

Conversion of Mass: 


3x ey ^ 


( 2 - 2 ) 


Conservation of Momentum: 


Du 
Dt " 

9u 

at 


3u , 3u „ 

“ toT + ^ - 

1 

Pf 

8p 

9x 


(2-3) 

Dv 
Dt " 

8v 

at 

+ 

9v , 9v 

u — + V -s — = - 

8x 9y 

1 

P£ 

9Pm 

By 

+ v^v + g Pj AT 

i 

(2-4) 


Conservation of Energy: 


Dt at ^ 9x ^ ^ 3y ^ ^ ^ c 


Zv . 


P 


P ^ (If)'] 
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VAPOR REGION ; 

Conservation of Mass: 


j 9 p V 


at 


+ 


8x 




9y 


= 0 


(2 


Conservation of Momentam: 


Du au , an , 0u 

5r = 8F + 




(2 


Dv 8v , 9v , 9v 
Dt ^ at ^ ax ^ By 


9Pr 


■t V V + gP^AT (2 


V 


Conservation of Energy: 


DT 1 °fPv^Pv> j. 1 

T 


Dt 


c Dt 

V 


P c 

V 


^^v , V^kVT 

T 


p c Dt P„ c 

r-y y f^V ^ 




(2 


Equation of State; 


= £(p^, T^) 

LIQUID- VAPOR INTERFACE BOUNDARY CONDITIONS : 
Velocity; ' 

Continuous Tangential Velocity: 





6 , 


- 6 ) 

-*7) 

- 8 ) 

-9) 

2 ” 

- 10 ) 

- 11 ) 


■ Normal Velocity; 


^gurface ^ S. * surface dt ^ j 


T emper atul e; 




Energy; 



o V, ^ 

Pf ^fg dt 


Shear; 



3 u ! 

9x 


dT dx 


EXTERNAL BOUNDARY CONDITIONS; 


Velocity: 



n(0,y, t) = v(0,y,t) = 

0 

; 0 S y S Y 

u(D,y,t) = v(D,y,t) = 

0 

; 0 < y <Y 

u(x, 0,t) > v(x, 0,t) = 

0 

; ' 0 < X £ L 

u(x,Y,t) = y(x,Y;t) = 

0 

; 0 X £ D 

Temperature: 


■ 


Xt ' 


9T I *“ Q 

ax " 1 0 (no ullage heating) 


(2-12) 

(2-13) 

(2-14) 

(2-15) 

(2-16) 

(2-17) 

(2-18) 

(2-19) 

( 2 - 20 ) 

( 2 - 21 ) 


1 


k |^(D, y, t) = q" : OSySL 


IJ(D.y.t) ={f : L<y^Y 


\ ^ •= - q” ; 0 S X <D 


k,|2{x,Y,t) ={g" : o<:.< 


D 


( 2 - 22 ) 
(2-23) 
(2^24) 
(2-25) 


INITIAL CONDITIONS; 


T(x„y, 0) = T^(x, y) 


u (x, y^ 0) = 0 


V (x, y, 0) = 0 


p{x, y, 0) = (x, y, 0) 


(2-26) 

(2-27) 

(2-28) 

(2-29) 


2. 3 CHARACTERISTIC VALUES 

■The characteristic quantities were chosen for the length, mass, temperatare, 
and time system and are; 

D - Characteristic Length (ft) 

pjgD3 - Characteristic Mass (Lbm) 

0 ‘ - Characteristic Time (sec) 

q"D/kjg - Characteristic Temperature (°R) 

q"0/D - Characteristic Pressure (p^) (Ib^/ft^) 


pje 


D2/q 2 „ Characteristic Dynamic Pre ssure (p^) {Ib^/ft ) 


where 

D is the tank diameter, 

■p^ is the density of the liquid phase, 
q” is the wall heat flux, and 

is the thermal conducti'wity of the liquid phase. 



The use of the tank diameter as a characteristic length and the liquid density 
to define the characteristic mass is a commonly accepted practice in model- 
ing, The selec tion of 0 as the characteristic, time is generally more conven- 
ient than VD/g. Here 0 could be cons.iderd^ ^p' be- the mission time or the 
stratification experiment time. The wall heat flux, q", multiplied by the 
ratio used to define the characteristic temperature, is also a value 

accurately established prior to most vehicle flights or experiments. The 
measured constant heater input heat flux, q", is used for determining the 
characteristic pressure and temperature values. 


2. 4 NONDIMENSIONALIZATION OF EQUATIONS 

The characteristic quantities are used to obtain the transformation of the 
differential equations from a dimensional to a dimensionless form. All 
variables can be nondimensionalized directly from these quantities except for 
pressure where the terms p^ and p^must be treated separately. The value of 

p is determined by fluid motion and should be nondimensionalized with a 

^2 2 2 
form such as pV . Hence, the value p^D /0 will be used. On the other hand 

p^ is essentially the measured tank pressure and is related to T through the 

equation of state. Thus, the selection of a characteristic value to nondimen- 

sionalize p^ can be obtained from the equation of state. It is convenient to 

use the perfect gas law 


p s p H.T 
V "^v V 


(2-30) 


which does not result in a loss of generality. Thus, nondimens ionalizing the 
right-hand side of the equation with the characteristic temperature yields: 


V ■ ra 


Pr R 

TZ “ 
P 


p T 

V 


(2-31) 


which indicates that the dimensional form (q'*0/D) can be used as a character- 
istic pressure, hence ’ 


Pv "" {q”0/D) 


V 


(2-32) 


g 


The form {q"9/D) can also be obtained through a physical argument where 
Pfinal (QipQip/m) «« (q”0/D). The resulting dimensional variables will be 


used 


X 


u 


X 

D 


u 


D/0 


m 


m 


T = 


p,D^/9^ 


±_ 

Q 




V = 


V 


AT = 


D/0 

Pv 


Wd 


AT 

q'«D/k, 


(2-33) 

(2-34) 

(2-35) 

(2-36) 


where 


9 = Time (duration of stratification test) 

NORMALIZED DIFFERENTIAL EQUATIONS: 
Liquid Region ; 

Conservation of Mass: 


9 u , d V 
9x 0y 


0 


(2-37) 


Conservation of Momentum: 
X component: 


Du 9u , — 8u , ~ 8u 

d7 = 87 + “ K + ^ 


35 


m 


9x 




111 


2 

V u 


(2-38) 


or 


Du 
Dt “ 



8x 


+ Fo'V^u 


(2-39) 



where 


Fo = 


(2-40) 


Y component: 


Dv 3v , — 8v . — 3v 
Dt 9t- 8x 9y 


.. AT 

^ V V+ 


9Prv, 7- . 2 — 

' + Fo • V V + Gr^ * Fo • AT 

By 


where 


Gr^' = 


gP^q" 

2 , 

'‘i h 


Conservation of Energy: 


DT _ ^ . •- BT ■ - 3T 
Dt “ Bt 9x ^ By 


V + -2 

d 2 q"° V 


' _ 2 , 2 

/b u\ ^B v j 


where 


so that 


= X£ 1 ' ^ . 2 fe? + fel 

Dt ” Pr Gr^ Fo c \dx./ \^y/ 
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Vapor Region : 

Conservation of Mass: 


3p 9pn 9pv 

V . . . . *-Y_ 

8t 9x By 


Conservation of Momentum: 


Du 9n , — Bu ,_9u 

Dt = sV “ iF ° 


, Bp V 0 , 

v^n 


„ 4- Fo-tx . v^u < 

ax 

S + u|x + v|5 = . V^v + ATI 

,8t BX By By 


^ Fo - V^v H- Fo^ P AT 
Pv Pv . "" 


Conservation of Energy: 


DT 2 _i 

■ Pi% D 


% Dt <Pv/Pv) >V p.c Pv 


\ y + /- 

c^q"D0 j \Bx/ \By/ J 


V + mY - -( ^ + ^Y 

X By/ 3 \ 8x 9y/ 


( 2 - 5 Z) 



or 


DT 
D T 


Fo 
Pr c 


V L 


D Y I 1 ^ 

Dt (Pv/Pv) + - £)t. 

ir'v 


+ |0 ^ ZjT 

Pr p c 

V 


(Z-53) 


■ -P-,r ■ 

1, 

faaY (sv'f 

M. Y§i j 

— Y 

2 /a u 

— Y 

Gr^ Fo n c 

t'v V 


\w) * \w) . 



“ 3 \3x ■ 



NORMALIZED LIQUID VAPOR INTERFACE BOUNDARY CONDITIONS : 
Velocity: 


TT. = tr 

j2 V 


{ 2 - 54 ) 


'^surface ~ 


V = V 


+ 


dy ~ 


V surface dr \ p 


(2-55) 


Temperature: 


T = T. 
V i 


(2-56) 


Energy: 


3 T 

9 y 






-IV 


qi >0 


|I dy 

dr 


(2-57} 


Shear: 


9u 

9 x, 


_ V - aa 

V jg v 8 x J 


q'^9 dcr dT „ Ma* dT 

p^v/^k^ dT dx ~ Pr dx 


( 2 - 58 ) 


NORMALIZED EXTERNAL BOUNDARY CONDITIONS : 
Velocity: 


a( 0 ,y,T) = v(d,yjT) = 0 ; Osy<Y/D 


(2-59) 



u(l,y,T) = 

v(l,y,T) = 0 

; 0 S y ? Y /D 


(2-60) 

u(x, 0,t) = 

v(x, 0, t) = 0 

; 0 < X < 1 


2-61 

u(x, Y/D,t) 

= v(x, Y/D, t) 

= 0 ; . 0 < X < 1 

> 

(2-62) 


T etnp e r atur e; 


|f(0,y,T) 

= - q”/q" 

= 1 

; 0 5 y < L 

(2-63) 


.. ( (- q"/q‘’ 
" |(0/q'') 

) = -1 
= 0 

L < y < Y /D 

(2-64) 

9T n - 1 

a=-(i,y,T) 

= q"/q” = 1 


Vi 

!!>' 

VI 

o 

(2-65) 


_ 1 (q'7q‘') 
"1 ( 0 /q”) 

= 1 
= 0 

; L < y < Y/D 

(2-66) 

|f(5,0.T) 

- q'7q'' = 

1 

; 0 < X S 1 

(2-67) 

^|f(5,Y/D,T) 

_ 1 (q'7q") 
- (0/q”) 

= 1 
^ 0 

; 0 < X < 1 

(2-68) 


NORMALIZED INITIAL CONDITIONS : 

T{x,y,0) = T.(x,y) (2-69) 

a{Xj y, 0) “0 (2-70) 

v{x,7, 0) = 0 (2-71) 

y 

p(xj y, 0) = Pj,(x, y, 0) (2-72) 
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2. 5 ‘ RESULTING DIMENSIONLESS PARAMETERS 

The dimensionless coefficients associated with the terms in the differential 
equations and boundary conditions are dimensionless parameters. Some of 
these coefficients include more than one of the dimensionless groups. When 
this is the case, the group is broken into more familiar subgroups, each of 
which are also dimensionless. The resulting set of dimensionless groups, 
along with their original forms, are presented in Table 2-1. 

Since it is not possible to scale two tanks while satisfying all the independent 

•JU aJU 

dimensionless groups, Gr*^, Fo, I, Ma"*, Pr and Cp/gp^D, two were elimi- 
nated. The first was the Marangoni number associated with motion at the 
liquid vapor interface resulting from surface tension changes with tempera- 
ture. Since little variation in temperature along the interface was expected, 
and because no such motion was ever observed during flow visualization 
experiments with Freon 113 conducted earlier at McDonnell Douglas, the 
Marangoni number was not kept. The other pi-group, c /g(3,D, was elimi- 

It 

nated since it is associated with the viscous dissipation term in the energy 
equation and is assumed to have a negligible effect during stratification. 
Table 2-2 gives the resulting set of dimensionless scaling groups that were 
used in this study. The dimensionless property ratios were all satisfied 
since the scaling tests were conducted with a single fluid. 
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Table 2- 1 




Remarks 

















Table 2-2 


SCALING GROUPS FOR THE COMPONENT, TWO-PHASE 
SYSTEM WITH A FLAT INTERFACE 


Parameter 


Requirement When Same Fluid 
Used for Model and Prototype 


2 , " 2 
D model D prototype 




gq”D' 


[model 


= gq"D^ 


Satisfied 


Satisfied 


Satisfied 


Satisfied 


Satisfied 


P^(t/8) 


Satisfied 


q"Dc, 


qD - qD 

model prototype 


Section 3 - 

EXPERIMENTAL EQUIPMENT AND INSTRUMENTATION 

3. 1 INTRODUCTION 

The stratification and de stratification experiments were conducted with a 
thermal system specifically designed for that purpose. The experimental 
system included three test tanks capable of supplying heat through the walls, 
a centrifuge facility for high-g tests, and the instrumentation necessary to 
regulate and monitor the experiment. 

The test tank included segmented blanket heaters covering the entire external 
surface of the tank and encapsulated by foam insulation. Thermocouples were 
used in the liquid, on the tank walls, on tank heat shorts (penetrations), and 
in the insulation to monitor the temperatures. In addition, a centrifugal 
pump was used for the destratification subsystem. The MDAC Santa Monica 
centrifuge facility was used for the high-g experiments and was necessary to 
achieve the same Grashof number for scaling the 6 inch diameter tank to the 
12 inch and 18 inch diameter tanks. The test liquid, Freon 113, was selected 
after evaluating a substantial number of candidate liquids. 

There were five major instrumentation components, 

1. Thermocouples for obtaining temperatures, 

2. A pressure transducer, 

3. A flowmeter for monitoring the destratification jet velocity. 

4. A variable power supply for each heater blanket. 

5. A centrifuge facility for producing the high-g levels required. 

Each of these items will be discussed further in the foUowing subsections, 

PEEGEDINa PA(JE BIAOTC NOT FILMED’ 
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3. 2 TEST TANK 


3. 2. 1 Basic Design 

The baseline Modular Nuclear Vehicle (MNV) and Saturn S-IVB-203 was 
selected as a typical tank and provided the general guideline for establishing 
tank geometry, tank dimensions, and the thermal properties associated with 
the test tanks. Hence, the test tank was a right circular cylinder with 
hemispherical dome caps on each end of the cylinder. De stratification sys- 
tem design is less well established so that the mixing system was designed 
as simply as possible using a system somewhat similar to that used in the 
■eariy basic de stratification experiments conducted at General Dynamics 
(2, 3). 

Three similar tanks of different dimensions were selected for the experiments 
since a single fluid was used. These three tanks in conjunction with the 
centrifuge facility were required for the preservation of the dimensionless 
parameters needed for similarity. These test tanks were fabricated from 
347 stainless steel with diameters of 6, 12, and 18 inches. The length of the 
tanks were 18, 36, and B4 inches, respectively, giving a length to diameter 
ratio of 3 for each. A detailed design drawing for the 6 in. diameter tank is 
shown in Figure 3-1. The dimensions shown were multiplied by a factor of 2 
and 3 for the 12 in. and 18 in. diameter tanks, respectively. The tanks were 
fabricated and designed to maintain geometric and thermodynamic similarity 
to satisfy the scaling parameters listed in Table 2-2, In addition to the 
scaling of the external dimensions, all wire diameters, wire insulation 
thicknesses, foam insulation, and wall thicknesses were scaled. 

It should be noted that the tank design using flanges at the junction of the 
hemispherical ends and the cylindrical section increases the versatility of 
the apparatus for further tests without compromising this study. In addition 
to allowing convenient access to the interior of the tanks to repair or relocate 
instrumentation, and to adjust the simulated jet mixing device, the tank end 
configurations could be altered by simply constructing new end sections and 
attaching them to the standard cylindrical section. The design also allows 
the ends to be connected to form a spherical test tank, with the cylindrical 
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section removed. All of these alterations could be accomplished while 
maintaining the geometric similarity between test tanks. Figure 3-2 is a 
photograph of these tanks. 

3, 2. 2 Test Tank Heaters 

An electrical heater blanket covering the entire tank was used for generating 
the various wall heat-flux values in the stratification experiments. The 
blanket made by embedding a resistance wire mesh into a thin silicone rubber 
sheet was produced by Electrofilm Inc. , North Hollywood, Calif. 

Ten heater blankets, each with separate electrical leads, were employed and 
had the capability to heat any part of the container, if desired. Thus, liquid 
heating, ullage heating, or a combination of the two could be achieved with 
this apparatus. Figure 3-3 shows the heater blanket arrangement; in this 
study the two heater sections at each elevation were wired in parallel to 
serve as one heater. This resTolted in five separate heater circuits. The 
heater blankets used on all three tanks were capable of producing a maximum 
of 10 watts/in. ^ (4, 910 Btu/hr-ft^). The blankets are shown in Figure 3-2. 

3.2.3 Test Tank Insulation 

Polyurethane foam was selected for the test-tank insulation because of its 
low thermal conductivity (0, 01 Btu/hr-ft°F), its relatively high maximum 
operating temperature {300°F), and because it could be foamed in place and 
is rigid after curing. The tanks and the heaters were essentially encapsulated 
within the foam. This was accomplished by placing the tank in a plywood box 
and filling the void between the tank and the box walls with foam. In this 
manner the plywood box became the protective covering, allowing for easy 
handling of the tank and providing a mounting platform for the tank during the 
stratification exp’.riments. In addition, smce the foam itself provided the 
structural support for the tank, heat short problems due to structural 
members attached directly to the tank were eliminated. Figure 3-4 shows 
the tank encapsulated by the foam in the plywood box. 

3. 2. 4 Pest ratification System 

The basic purpose of the destratification system was to mechanically and 
thermally mix the cold liquid near the tank bottom with the warmer liquid 
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FIGURE 3-4 

FOAM ENCAPSULATED TEST TANK 



above it to study its effect on the tank pressure. An axially directed nozzle 
located in the bottom of the tank and connected to the outlet of the centrifugal 
pump was used to circulate the cold liquid up through the warmer liquid 
region. An inlets concentric to the nozzle and leading to the pump inlet, 
provided for the closed-fluid loop that resulted when the pump was operating. 
In this manner, the cooler liquid along the bottom was fed into the pump and 
exited as a cold jet of liquid moving upward into the warmer liquid. 

Figure 3-5 is a sketch of this flow loop. 

The nozzle openings for each of the tanks were scaled in the same manner as 
the tanks. The nozzle openings for the 6, 12, and 18 inch diameter tanks are 
0.0625, 0. 125, and 0. 1875 inches in diameter, respectively. The inlet and 
outlet tubes protruding from the tanlc were also scaled. Exact scaling was 
not maintained with the external piping which connected the pump to the tank. 
However, the size and length of this plumbing increased approximately in 
proportion with increasing tanlt size. All of the external piping was wrapped 
in fiber glass insulation to reduce the heat loss to the surroundings as the 
heated liquid was circulated through the pump. An Eastern 1/15-hp centri- 
fugal pump was used for all three tanks. The flow rate for this pump was 
regulated with a Superior Electric Co. , 7.5 amp. Power stat. 

3. 3 INSTRUMENTATION 

3. 3, 1 Thermocouples 

z\ large number of thermocouples were necessai’y for the stratification study 
to accurately assess all of the transient thermal changes that occur during 
the heating and de stratification periods. Thermocouples were used to 
measure the temperature distribution within the liquid and vapor, and to 
monitor the heat stored in and/or lost through the tank walls, insulation, and 
the heat shorts. The wall and insulation measurements were needed to 
determine the energy stored and hence to provide a reasonably accurate 
estimate of the actual heat flux to the liquid in the tank. However, the num- 
ber of thermocouples was limited by the instrumentation, and by the risk 
of the thermocouple rake interference with the convection currents produced 
by the de stratification pump. A flow visualization study with a plexiglass 
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FIGURE 3-5 

DESTRATIFICATIOIS! FLOW LOOP 
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DESTRATPUMP 




tank and several thermocouple rake configurations was conducted to find 
the rake which had the least influence on the destratification pump flow 
patterns. 

Figure 3-6 is a schematic of the rake selected and the thermocouple 
locations. The two cross members allow for the radial temperature distri- 
bution at two levels. These rakes are also shown in the photograph of 
Figure 3-4. A total of 37 thermocouples are on the rake. The liquid and ull- 
age thermocouple locations relative to the tank height are given in Table 3-1, 
Figure 3-7 shows the thermocouple attachment points on the tank walls as 
well as on penetrations and within the foam insulation. A total of 40 thermo- 
couples were used to monitor the wall and penetration temperatures. 

The thermocouples were made from gauge number 36, 30, and Z6 copper 
constantan wire for the 6, 12, and 18 inch diameter tanks. Two, 50-terminal 
Research, Incorporated, 150 “F thermocouple reference junction ovens 
were used with the thermocouples. The output of the 77 thermocouples was 
handled with a Hewlett Packard EOlOB Data Acquisition System. A 291 IB 
Hewlett Packard Crossbar Scanner Control was used to scan and print out 
all 77 voltage values in less than 15 seconds. The thermocouple signal was 
amplified with a Dymec Guarded Data Amplifier and then sent to the Hewlett 
Packard 240 1-C Integrating Digital Voltmeter. The voltages were printed 
with a Hewlett Packard Model J66-562A Digital Recorder. The scanner was 
linked to a Hewlett Packard digital clock so that a scan frequency could be 
determined and set prior to the test period. This data acquisition system is 
shown in Figure 3-8. A schematic for the entire test arrangement is shown 
in Figure 3-9. 

3.3.2 Tank Pressure 

The tank internal pressure was monitored and recorded along with the thermo- 
couple data. A pressure transducer, located at the top of the tank, also had 
its output signals recorded in the Data Acquisition System. A Statham 0 to 
50 psia, 7-volt, maximum-pressure transducer was used for the autoxnatic 
monitoring of the pressure. A Trygon Electronics Model HR-20- 16A Power 
Supply was used with this transducer. Observation of the pressure history 
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Table 3-1 

LIQUID-ULLAGE THERMOCOUPLE - NORMALIZED TANK HEIGHT RATIOS 
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with a visual pressure gage was made along with the transducer recordings. 
The visual pressure gage used for the IG tests v/as a Seegers Precision 
Pressure Transfer Standard, marked in 0. 10 psia subdivisions. 

3.3.3 Blanket Heater Power 

The blanket heaters were simple electrical resistance heaters. The five 
heater pairs (1-2, 3-4, 5-6, 7-8, 9-10) were operated separately, each 
using a General Radio Type WZOMTS 20 amp Variac to regulate the voltage 
in a continuous manner. The power for each heater segment was measured 
with a Singer Single-Phase, Low-Power Factor Wattmeter, model DLW, with 
a range from 150 to 1,200 watts. For the experiments where the power 
exceeded 1, 200 watts, a separate ammeter and voltmeter were used to 
measure the power output. 

3. 3. 4 Flow Meter for Destratification Pump 

The mass flow rate of the de stratification jet was monitored by measuring 
the volumetric flow rate of the F reon on the outlet side of the destratification 
pump. A Waugh Turbine Flow Meter (Instrument No. 101877-102-1327) with a 
flow range from 0. 15 to 1. 5 gpm was used. The output of the flowmeter was 
fed into a Waugh Pulse Rate Converter, Model RF-213. Frequency values 
from the converter were then translated into flow rate values. 

3.4 SIMULANT LIQUID SELECTION 

The first stage of the screening process considered 40 noncryogenic liquids 
and Reference 4 documents that process in detail. The criteria considered 
to be important for scaling tests with a single fluid were: 

A, Boiling point of the liquid. 

B. Ease of handling. 

C, Availability of precise property data. 

D. Cryogenic propellant similarity. 

The boiling point of the liquid, was an important consideration since a liquid 
which boils near room temperature would cause fewer problems in general. 
For example, a low boiling point results in less heat loss to the surroundings, 
make handling easier (as well as test control), and also makes it easier to 
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reach an isothermal condition. In addition, it is convenient to have the room 
temperature and saturation temperatures close together. This reduces the 
time between experiments whenever the tank must be heated or cooled down 
to the saturation temperature. The availability of accurate property data 
was mandatory to interpret the experimental results. 

Table 3-2 lisiis some of the readily available liquids with boiling points 
between 94, 3°F and E31°F, Flammability and toxicity were also considered 
to be important factors along with availability. As a result. Freon 113 PCA 
(i. e. , Precision Cleaning Agent) was selected as the liquid best suited for 
this program. The boiling point of Freon 113 is approximately 117°F and as 
a result does not cause burns when it contacts the skin. Freon 113 is also a 
relatively safe liquid to use with proper ventilation. Threshold limit values 
(TLV) adopted by the American Conference of Governmental Industrial 
Hygenists provided a guide to the toxicity of solvents and other substances. 
These limits in ppm of air for Freon are the same as for ethyl alcohol and 
acetone and are at the highest TLV for any material with the exception of 
carbon dioxide. In addition, the fluorocarbon solvent vapors are completely 
nonflammable, Th? Freon 113, also labeled as Precision Cleaning Agent 
(PCA), is relatively pure ther.-* ^ minimizing solid deposition within the 
tank since makeup liquid was required after each test. 

3. 5 CENTRIFUGE FACILITY 

The MDAC-W 3 6-ft- diameter centrifuge facility was used in this study to 
obtain the 8-g and 27-g environments for the thermal stratification tests. 

The large centrifuge was used to minimize the gravity gradient across the 
6 inch tank. This is the result of the fact that the acceleration due to rotation 
is ro}^. Because of the 18 ft arm and the 18 inch tank length, the change in 
the g level along the tank axis is approximately 10 percent. The tests were 
performed with the 6 inch tank using the same instrumentation that was used 
on the 12 and 18 inch tanks 1-g environment tests. 

The centrifuge was equipped with slip rings to permit application of electric 
power for heating the test tanks and monitoring the thermal effects in the 
test fluid. Figure 3-10 shows the tank mounted on the centrifuge. 
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Table 3-Z 

CANDIDATE FLUID BOILING POINTS 


Fluid 

Boiling Point (°F) 

Acetone 

133.2 

Ether (diethyl -) 

94. 3 

Pentane 

96.8 

TrichlorotriUuoroethane 

117. 6 

Chloroform 

142.2 

Water 

212.0 

Hexane 

154, 5 

Ethyl Acetate 

77. 1 

Toluene 

231. 1 

Methyl Alcohol 

149.0 

Benzene 

176.2 

Heptane 

209. 1 

Carbon Tetrachloride 

170. 2 

Ethyl Alcohol 

173.3 

Allyl Alcohol 

206, 6 


The centrifuge is powered by a closed-loop hydraulic drive system of 
450 horsepower which urovides control to within 1/10 rpm. The centrifuge 
arms, measured from the center of rotation/ are of unequal length. Both 
ends of the centrifuge are commissioned to accept fixtures and instrumenta- 
tion wiring. 
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Various configurations of closed- circuit television are available for test 
monitoring. Operational capabilities of the centrifuge are listed below: 


Manufacturer Model 

Rucker 58-3189 

Rated capacity 

76, 000 1b 

Maximum accelerations 

150 g’s 

Maximum specimen weight 

30 0 lb 

Pit diameter 

41 ft 

Long-arm radius 

15 ft, 10 in. 

Short- arm radius 

11 ft, 4 in. (18 ft with 
platform extension) 

Height of arm centerline 
above floor 

4 ft, 4 in. 

Slip rings (175) Power 

4 at 200 amps 
9 at 20 amps 

Instrumentation 

16 at 6 amps 
146 at 2 amps 

RF 

1 coaxial connector 

Universal joints 
hydraulic or 
pneumatic 

5, 000 psi (2 systems) 
2, 500 psi (2 systems) 
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Section 4 

TEST PROCEDURE 

4. 1 INTRODUCTION 

The experimental procedure for the stratification experiments consisted of 
three major steps: 

1. Conditioning the tank, liquid, and vapor. 

2. Thermal stratification test. 

3. De stratification test. 

The conditioning procedure included the steps necessary to reach the 
saturation condition for the Freon 113 and its vapor at atmospheric pressure. 
In addition, it was necessary to ensure that any liquid motion produced by the 
mixer was damped out prior to the initiation of the stratification test period. 

t 

The thermal stratification period was established by practical limits on the 
maximum allowable tank pressure and the maximum allowable temperature 
of the electric wall heaters. This period also had a lower bound associated 
with the time to reach a quasi- steady- state temperature profile. The 
destratification period was established from earlier destratification studies 
(Ref 5). These results were used to define the de stratification test period 
for this program. Details for performing these major steps are included in 
this section. 

4. Z LIQUID, VAPOR, AND TANK CONDITIONING 

Conditioning of the Freon 113 prior to the start of a stratification test was 
necessary to meet the following conditions: 

A. An isothermal system at the Freon 113 saturation state (i. e., at 
one atmosphere), 

B. The liquid and vapor must be quiescent. 

The system was brought to the saturation state with the aid of the wall 
heaters. The tank vent valve was kept open so that the Freon reached the 
saturation temperature at one atmosphere. The de stratification mixer was 
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left on during this heating period to aid in maintaining the tank walls and 
liquid at a uniform temperature. At this point, liquid motion due to the mixer 
and convection was present. This residual motion was then allowed to decay 
prior to initiating the stratification tests. 

The required decay period was determined in an approximate fashion by 
operating the de stratification pump in the 12 inch tank with the upper dome 
removed. The pump was turned off and the time required for the visually 
detectable residual motion to cease was noted. This time was of the order 
of two minutes and this value was used as a rough guide for the delay period 
prior to the initiation of a stratification test. 

The same steps were taken during the centrifuge experiments. Here, the 
conditioning operations were conducted with the tank in the vertical position. 
This was necessary for normal heating and venting, where just before the 
centrifuge was put in motion the test container was positioned on its side so 
that its axis remained parallel to the resultant gravity vector of the designated 
''g" level. 

4. 3 STRATIFICATION TEST PERIOD 

Several criteria were considered when determining the thermal stratification 
period, 0 . One criterion was the time required to reach a quasi- steady- state 
temperature profile in the liquid. The other stratification test period 
estimate was the time required for the mean temperature of the heated liquid 
to rise 15°F, 

Chin et al (Ref 7), Neff (Ref 8), Schwartz and Holmes (Ref 9), and Tellep and 
Harper (Ref 10) have all reported that stratification profiles appear to main- 
tain the same shape after an initial transient period. As a consequence, when 
a quasi- steady- state condition is reached, the propellant tank should act as a 
lumped-mass system and, fchepressure rise rate can be reliably predicted. 
The time period using this criterion was obtained from a semi- empirical 
equation obtained by Holmes (Ref 6) using Freon 113. Although the tank 
used by Holmes did not have the same geometry as the tanks in this study, his 
equation was selected for comparison with the 15°F temperature rise criterion 
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as a first approximation. Holmes (Ref 6) obtained the following correlations 
to fit his data: 

Fo^Ra^^^^ = 0.2“ (laminar) (4-1) 

and 

Fo»Ra^^^*^ = 1,86“ (turbulent) (4-2) 

14 1 6 

The values of Gr^j in this study ranged from 10 to 10 which is generally 
considered turbulent. However, the lower portion of the tank walls will 
still be laminar. Therefore, to establish the longest test time required with 
this criterion the laminar equation was used. The overall system balance 
was obtained by equating the heat added to the liquid over the stratification 
test period to the sensible heat gain of the bulk liquid. Here 

q"A^B =Q=xn^Cp^(T^^.T^,_) (4-3) 

where T-j^^ - Tf^^ ~ 15‘^F and A.^, is the wetted area of the tank. It was found 
that the values of 0 obtained from equation (4-1) were nearly the same as 
those obtained using the 15°F mean temperature rise requirement. As a 
result, the 15°F rise criterion was used to determine 9. 

4. 4 DESTRA.TIFI CATION TEST PERIOD 

The de stratification test period was initiated by actuating the simulated jet 
mixer. Since it was of interest to operate the mixer with and without the 
wall heaters on, both configurations were considered for some of the experi- 
ments. When the wall heaters were left on, the mixer had to overcome the 
free convection currents in addition to the buoyancy of the stratified fluid. 

The thermal de stratification restolts of Poth (Ref 3), without wall heaters 
operating, clearly indicate that axial jets are superior to radial jets; there- 
fore, only axial jets were used in this study. The results of Poth have also 
been used to define the minimum destratification time for the unheated cases. 
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The upper bound of data presented by Poth recommends the following 
conservative estimate of the mixing time period: 


® M “ V.D.’ 

J J 


where 


e 


M " 
D ' 


V. = 
J 

D. = 
J 


Mixing time 

Tank diameter 
Jet velocity 

Jet diameter 


(4-4) 


It should be noted that the predicted values were only considered as a guide 
since monitoring of the temperature data during the testing provided the 
most reliable assessment of when the de stratification period should be ter- 
minated. The time required for mixing with the heaters in operation was not 
determined prior to the test. Hence, the mixer was kept in operation until the 
temperatures in the liquid and the tank pressure indicated that mixing had been 
achieved or until such time that it became apparent that the jet mixer would 
not overcome the prevailing convective fluid currents. 


An analysis was conducted to provide a reasonable upper limit on the jet 
Reynolds number, A critical jet Reynolds number was defined as that value 
where the buoyant force becomes unimportant when compared to the inertial 
forces. For a Reynolds number greater than this critical value, the system 
may be assumed to mix completely and increasing the jet Reynolds number 
further simply decreases the mixing time required. Below this critical Rey- 
nolds number, the buoyant forces may be strong enough to limit the degree 
of mixing. 


Although a number of analyses were performed to determine this critical 
Reynolds number, the maximum value is given by 



(4-5) 
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where 


Re. - u.d/v-^ ^ 

J J 

u. = Jet velocity at orifice, and 
J 

d = Orifice diameter. 

This criterion was established by assuming that the mixing jet must possess 
enough energy at the free surface to overcome the buoyant force, and move 
the hot liquid at the surface to the bottom of the tank. Using this jet Reynolds 
number as the critical value, the mixing jet should easily penetrate the strati- 
fied layer where the momentum due to the pump is always greater than the 
free convection momentum, and the liquid kinetic energy due to the pump is 
always greater than the free convection kinetic energy. 

This latter equation was only used as a guide to place a reasonable upper 
limit on the jet Reynolds number requirements for this study. Observations 
of the temperature data during de stratification also were used to ensure that 
the critical number was reached. 
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Section 5 

EXPERIMENTAL PROGRAM DEFINITION AND TEST MATRIX 
5. 1 STRATIFICATION 

The test matrix presented in this section -was defined to check the validity of 
the dimensionless parameters chosen in Subsection Z. 5. Here, scaled pairs 
of experiments were selected for various combinations of the three different 
sized tanks. These pairs were arranged to give the widest possible range of 
values for the dimensionless numbers. Thus, if the dimensionless pressure 
and temperature values for each scaled pair of tests match each other then 
these dimensionless groups would be adequate to scale the stratification 
process. Limitations on the range of values for the dimensionless groups 
are discussed below. 

The first step in defining the bounds on the experimental program was to 
limit it to that of a single component. Freon 113, and its vapor. Three 
scaling parameters, the Modified Grashof number (Gr ), the Fourier num- 
ber (Fo), and the Interface number (I), were kept constant for the tests which 
were scaled one to another. Scaling tests were conducted between the 6 inch 
diameter tank and either the IE inch or 18 inch diameter tanks. The experi- 
ments selected are listed in Table 5-1, All experiments were conducted with 
liquid heating and some with simultaneous ullage and liquid heating. The test 
matrix was designed to cover a range of Modified Grashof numbers from 10^^ 
to 10^^. These limits were established by the Eillowable heating rates. The 
upper limit on Gr"^ was set by one of the assumptions made in the dimensional 
analysis. This was necessary since the governing equations did not include 
the effects of nucleate boiling on the tank walls. High heat-flux values in 
the vigorous nucleate boiling range could invalidate the descriptive 
equations and the resulting dimensionless parameters. For heat flux 
values in excess of 1, 000 Btu/hr-ft , boiling is expected at atmospheric 
pressure. As pressurization occurs, boiling should be somewhat depressed; 
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Table 5-1 

TEST MATRIX - FREON 113 {PCA} 


Test Column H 1 

Test Column tfZ 



Test Column #3 

6 in. Dia Tank 

12 in. Dia Tank 

18 in, Dia Tank 

aC - Acceleration - 27G 

IG 

1C 

1,000 Btu/hr ft^ (L) 87.5% FIJ_L 
Cr*=lCll5 TESTS #Z,i?l2S 


500 Btu/hr ft^ (L) 87, 5% FILL 
Gr=»= i0l5 TESTS #30. #25 


1. 000 Btu/hr tL+U) 87. S% FILE 

Gr* = IQ15 tests «3, #115 


500 Btu/hr (L+O) 87. 5% FILL 

Gr'i'3 10l5 tests #31. #17 


700 Btu/hr it^ (L) 50% FILL 

Grs= = 10^4 tests jgs, ffiss 


350 Btu/hr It^ (L) 50% FILL 
Cr«= = 10 14 TESTS #33, #32 


700 Btu/hr ft^ {LtU) 50% FILL 
Crf = 10 14 TESTS fl lOS, ffl6S 


350 Btu/hr ft^ (L+U) 50% FILL 
Gr=f = 10 14 tests #34, #35 


600 Btu/hr ft^ (L) 87. 5% FILL 
Gr=? = 7 3c l0 14 TESTS 111, 1/5, i/BS 


300 Btu/hr tt^ lU 87. 5% FILL 
Gr¥ = 7 X 10l4 TESTS 920, #22 


600 Btu/hr ft^ <L+U) 87. 5% FILL 
Gr^- = 7 X 10 14 TESTS ff4, #6 


300 Btu/hr ft’ (L+U) 87. 5% FILL 
Gr« = 7 X 10l4 tests #23, #28 


90 Btu/hr ft^ (L) 87. S% FILL 
Gr* = 10 14 TEST #14 


45 Btu/hr ft^ (L) 87, 5% FILL 
Gr* = 10l4 TESTS #26, #15 


90 Btu/hr Ct^ (L+U) 87. 5% FILL 
Or=."- I 10 14 test #13 


45 Btu/hr ft^ (L+U) 87. 5% FILL 
Gr=)== 10 14 test #27 



2,230 Btu/hr ft^ (Lj 50% FILL 
Gr* = 10l5 test #1S 


748 Btu/hr ft^ (L) 50% FILL 
Cr=:=^l0l5 tests #BS,#D 


2, BBO Btu/hr ft^ (L) 87. 5% FILL 
Gr'f = 10 lo TEST #3S 


960 Btu/hr ft^ (L) 87. 5% FILL 
Gr* = 10 16 TEST #1S 


NOTES: 

• (i.4U] « . . Liquid, and Liquid 4 Ullage Heating. 


• Testa horiaentally acroas Table are scaUng tceta ii*c. ^ tests having conatant Modified Grashof no. Gr* - = f (L = liquid depth)). 

• "S*' ... S in tests number indicates no heating during destratification testing. 
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however, a quantitative assessment of this effect is not known. Nevertheless, 
several experiments were conducted at heat-flux values near 3, 000 Btu/hr-ft , 

This value of heat flux was selected as a conservative upper limit for protect- 
ing the heater blankets. Therefore, it was not anticipated that the effect of 
boiling on scaling would be specifically determined other than by looking for 
changes in the pressure rise rate and the stratification profiles in the 
higher heat-flux range. The lower bound on Gr ^ or in effect q", was the 
degree to which the tank insulation acts as an adiabatic wall. Rough 
estimates on the heat leak through the insulation with the 12 inch tank were 
made from preliminary data and indicated that this value was between 1 and 
10 Btu/hr-ft^. Thus, 45 Btu/hr-ft^ was arbitrarily selected as the lower 
limit on q”. This condition corresponded to a Modified Grashof number of 
10^4 and it was the lowest value studied. 

Specific values of the Fourier number would not be an important consideration 
other than it be in keeping with expected values for space vehicles. Thus, to 
check scaling, a convenient value of the Fourier number was selected. Here, 
the stratification period 9 was selected to obtain an increase in the liquid 
temperature of the order of 15*^F which also corresponded to a pressure rise 
in the tank of the order of 5 psig. Smaller values could result in serious 
experimental error since temperature increases should be approximately 
an order of mangitude greater than the thermocouple error and larger values 
could tend to reach the safety limits of the tanks. The third dimensionless 
parameter did not include any new terms associated with the tank and thus 
did not impose any further limits on the experimental range of variables. 

The fraction of the tank volume occupied by the liquid in this study had two 
values, 50 percent and 87. 5 percent. Scaling between two tanks required 
that the percent fill be the same. 

5.2 DESTRATIFICATION 

Although the scaling procedure did not specifically include the de stratification 
process, there is reason to believe that scaling during this period should 
require the same dimensionless numbers as used in the stratification period 
along with the Reynolds number. As a result, the jet Reynolds number was 
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treated as an additional pi-group during this period. Table 5-E lists these 
jet Reynolds numbers; and, it should be noted that all scaled tests were not 
carried over in the destratification test region. 


Table 5-2 


TEST MATRIX - DESTRATIFICATION JET REYNOLDS NUMBERS 


Test Column #1 
6 in. Dia Tank 
8G - Acceleration - E7G 

Test Column #2 
12 in. Dia Tank 
IG 

Test Column #3 
18 in. Dia Tank 
IG 

Tests #2, #12S 

Re. s 38K, 38K 
J 


Tests #30, #25 
Re. = 38K, 3 IK 
3 


Tests #3, #liS 

Re. = 38K, 38K 
J 


Tests #31, #17 
Re. = 38K, 3 IK 
3 


Tests #9S, #15S 

Re. = 38K, 38K 
3 


Tests #33, #32 
Re. = 3 IK, 24K 
3 


Tests #10S, 16S 
Re. = 38K, 45. IK 


Tests #34, #35 

Re. = 3 IK, 20K 
J 


Tests #1, #8S, #5 

Re’. = 3 IK, 38K, 31K 
J 


Tests #20, #22 

Re. = 3 IK, 3 IK 
3 


Tests ^t4, #6 

Re. = 3 IK, 3 IK 
3 


Tests #23, #28 
Re. = 3 IK, 3 IK 


Test #14 

Re. = 38K 
3 


Tests #26, #15 
Re. = 3 IK, 25K 
J 


Test #13 

Re. r 38K 
3 


Test #27 
Re. = 31K 
3 


• 

Test #1S 

Re. = 38K 
3 


Test#BS, #D 
Re. = 38K, 38K 
3 


Test #3S 

Re. = 45. IK 
3 


Test#IS 
Re. = 45. IK 
3 


Note: Re. . 

J 


Reynolds number of de stratification (mi^dng) flow at 
nozzle exit. 
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Section 6 

DISCUSSION AND ANALYSIS OF RESULTS 
6.1 INTRODUCTION 

The data in this study correspond to those experiments listed in the test 
matrix (Table 5-1) for the 6, IE, and 18 inch tanks. Reproducibility is the 
first item covered to establish data reliability. Here, curves are presented 
for sets of identical experiments to indicate the spread in data due to 
experimental error. The Modified Grashof number definition was expanded 
to differentiate between the heat flux produced by the blanket heaters and 
that which actually enters the tank fluid. Then, the criteria which must be 
met to obtain scaling are discussed. Finally, Modified Grashof number, 
Fourier number, Interface number, pressure, and temperature data are 
presented and examined with these criteria to assess the degree to which 
scaling was achieved. 

The test liquid (Freon 113) properties, characteristic pressure and tempera- 
ture values, the values of the three scaling parameters; Modified Grashof 
numbers (Gr;^ and Gr^), Fourier number (Fo), and Interface number (I), for 
each test investigated £tre included in Appendix A. A summary of the heat 
absorbed by the liquid, the tank walls and heaters, stratification test time, 
and the heat source rates are tabulated for each test investigated in 
Appendix B. These tables show that the heat accounted for through thermo- 
couple measurements for all tests was about 95.5 percent of the h' "t actually 
generated by the heater blankets. This provided a measure with wiiich to gage 
the accuracy of the temperature values in the liquid, the ullage, and the tank. 

6. E DATA REPRODUCIBILITY " ) 

O.E.l Introduction 

Data reproducibility was evaluated according to the measured-bulk liquid 
temperature and tank- pres sure percent deviations achieved from repeated 
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test conditions conducted over a period of time ranging from days to months. 
Fourteen tests were repeated for reproducibility. The degree of reproducibility 
is discussed below. 

6.Z.Z Pressure History Curves 

The pressure history curves shown in Figures 6-1 through 6-14 correspond 
to some of the test listed in the test matrix {Table 5-1). These curves show 
good agreement for all test operating conditions in the Modified Grashof num- 
ber range of 10^^ through 10^^, for the range of 50 and 87. 5 percent liquid 
fill, and for the liquid, and liquid and ullage heating conditions. It was noted 
that the initia-l pressure in some tests were not the same. This is partially 
due to the variation in the barometric pressure conditions. However, in these 
test cases with the initial shift in pressure, the reproducible pressure 
curves have remained essentially parallel. Hence, the same increase in 
pressure during the repeated tests indicates that the variation of the pressure 
data from test to test was small. 


In addition to the graphical representation, a quantitative assessment of the 
pressure variatn’ns was made by determining the percent deviation of the 
re’- roducible pressure data. This was done by determining the pressure 
increase of each of the respective tests and then dividing it by the averaged 
initial pressure of each of the two reproducible tes-ts to be compared. 


This percent deviation is defined as follows: 

^PA-B _ PAf - Pa, - PBj^ 

Piavg " 1/2 (PA; + PBi> ■ 1/2 (PAi + PBi> ‘ 


The maximum deviation of pressure reproducibility is about 9.8 percent. 
Table 6-1 gives the values for all of the tests. 

Based on the reproducibility of the pressure data curves. Figure 6-1 through 
6-14, and the pressure deviation percentages shown in Table 6-1, it is 
concluded that a high degree of reproducibility was achieved between two or 
more identical tests. It is noted that these tests were repeated in a random 
fashion. Tests were repeated from 1 to 30 days after the initial tests were 
conducted. This further demonstrated not only the reproducibility and 
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PRESSURE REPRODUCIBILITY CURVES 
FIGURES 6-1 THROUGH 6-14 
NOMENCLATURE 

= Measured heater heat flux iuput, Btu/hr ft 
= Reynolds number at mixer jet nozzle 



TANK-ULLAGE PRESSURE tPSIA) 


CR 138-1 


FIGURE 6-1 

REPRODUCIBILITY DATA 

6 IN. TANK - PREiS. HISTORY: TESTS 8G NO. 2 AND NO. 12S 
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TANK-ULLAGE PRESSURE (PSIA) 
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CR 108-1 


FIGURE 6-3 

REPRODUCIBILITY DATA 

6 IN. TANK-PRESS. HISTORY: TESTS 8G NO. 9S AND NO. 15S 


q » 700 BTU/HR.-FT‘ (LIQ. HTG.) 50 % FILL 
H 
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Tank-ullage pressure (psia) 


FIGURE 6-4 

REPRODUCIBILITY DATA 

6 IN. TANK-PRESS. HISTORY: TESTS 8G NO. 10 AND NO. 16S 


CR 108-i 





FIGURE 6-5 

REPRODUCIBILITY DATA 

6 IN. TANK-PRESS. HISTORY: TESTS 8G NO. 1, NO. 5, AND NO. 8S 


CR 108-1 
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CR 108-1 


FIGURE 6-8 

REPRODUCIBILITY DATA 

12 IN. TANK-PRESS. HISTORY: TESTS IG NO. 31 AND NO. 17 
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CR 100-1 

FIGURE 6-9 

REPRODUCIBILITY DATA 

12 IN. TANK-PRESS. HISTORY: TESTS 1G NO. 33 AND NO. 32 
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CR108 I 


FIGURE 6-10 

REPRODUCIBILITY DATA 

12 IN. TANK-PRESS. HISTORY: TESTS 1G NO. 34 AND NO. 35 
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FIGURE 6-11 

REPRODUCIBILITY DATA 

12 IN. TANK-PRESS. HISTORY: TESTS IG NO. 20 AND NO. 22 


CR 1081 
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LLAGE PRESSURE (PSIA) 


CR108-I 

FIGURE 6-12 

REPRODUCIBILITY DATA 

12 IN. TANK-PRESS. HISTORY: TESTS 1G NO. 23 AND NO. 28 
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TANK-ULLAGE PRESSURE (PSIA) 


CR 108-1 


FIGURE 6-13 

REPRODUCIBILITY DATA 

12 IN. TANK-PRESS. HISTORY: TESTS 1G NO. 26 AND NO. 15 
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TANK-ULLAGE PRESSURE (PSIA) 











Table 6-1 


PRESSURE DEVIATION (PERCENT) 
(REPRODUCIBILITY) 



G?h 

Transducer. 

Pressure 

(Percent) 

Gage 
Pressure 
( Percent) 

6-Inch-Diameter Tank 





S0% Fill 

Liquid Heating 

10l4 

1.7 




Liquid-Ullage Heating 

10l4 

2.0 



87.5% Fill 

Liquid Heating 

7 X 10l4 

3.5 





10l5 

1.0 




Liquid-Ullage Heating 

7 X 10l4 

2.8 





10l5 

4.3 


12- 

[nch-Diameter Tank 




1 

50% Fill 

Liquid Heating 

10l4 

2.2 

1.7 

■ 


Liquid-Ullage Heating 

10l4 

6.0 

5.4 

■ 

87.5% Fill 

Liquid Heating 

1014 

8.0 

9.8 

1 



7 X 10l4 

1.0 

3.9 




10l5 

0.7 

1.3 

1 


Liquid-Ullage Heating 

7 X 10l4 

2.4 

2.8 

■ 



10l5 

6.7 

8.4 

18- 

[nch-Diameter Tank 




■ 

50% Fill 

Liquid Heating 

10l5 

1.0 

1.0 
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reliability but also the consistency of the entire test tank system models and 
^ apparatus over the range from 1 to 27 g's. 

6. 2. 3 Temperature Profile Curves 

The dimensional and dimensionless temperature profile curves at the 
normalized time (t) equal to 1.0 are shown in Figures 6-15 through 6-31. 

The curves show that there is good reliability in achieving reproducibility 
of the temperature measurements throughout the fluid. In addition, since the 
dimensionless temperature (T) and the dimensional temperature (T) only 
differ by a constant, both scales can be used to represent the same curve on 
any given graph. An examination of these temperature profiles for identical 
test conditions shows that the reproducibility tests nearly collapse on top of 
one another, with little difference in the temperature readings at each thermo- 
couple location. Here, the temperature scatter is less than 2°F at each 
thermocouple location. The reliability of reproducing test data over an 
elapsed period of time of 5 and 17 months apart is demonstrated by inspection 
of the temperature profiles shown in Figures 6-29 through 6-31. These 
curves show good reproducibility. 

The normalized liquid-bulk temperatures (average liquid temperature) shown 
in "Figures 6-32 through 6-45, were also analyzed for determining the percent 
deviation values (Table 6-2). The method used is identical to that used for 
determining the percent deviation values for the pressure data discussed in 
Subsection 6. 2, 2. 


Here, the percent deviation is defined and expressed as 

ATa_B ^-^f " ^-^i ^®f ' ^®i 

^Tiavg " ^Ai + TBi) ' 1/2 (TAi + Tb^) 


6-2 


Note that this technique automatically adjusts the comparison to account for 
the differences in the starting temperature whereas. Figures 6-32 through 6-45 
do not. As a result, when the starting temperatures are different, the curves 
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TEMPERATURE PROFILE REPRODUCIBILITY CURVES 


FIGURES 6-15 THROUGH 6-31 
NOMENCLATURE 

Measured heater heat-flux input, Btu/hr-ft^ 

”s" Symbol in test number indicates no 

destratification heating 

(L) Liquid heating only 

(L+U) Liquid and ullage heating 

Liquid Freon 113 (PC A) 

T = (t/0) where t = Normalized time 

t = Elapsed time test 
0 = Stratification test time 
6-In. Tank. . . . Title refers to tank diameter 

If 

Grpj Modified Grashof number based upon q^^ 
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CR 108-1 


FIGURE 6-15 

REPRODUCIBILITY DATA 
6 IN. TANK - STRAT TEMP PROFILES ( T = 1.0) 
TESTS 8G NO. 2 AND NO. 12S (L) HEATING 




CR 108-1 


FIGURE 6-16 

REPRODUCIBILITY DATA 
6 IN. TANK - STRAT TEMP PROFILES ( r = 1.0) 
TESTS NO. 3 AND NO. IIS (L + U) HEATING 



0 10 20 30 40 

DEL TMP = T (r) - T (t = 0) (°F) 


2 4 8 12 16 20 24 28 32 

NORMALIZED DEL TMP X 10^ 
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TESTS 8G NO. 9S AND NO. 15S (L) HEATING 








rANK HEIGHT 




0.4 0.5 


NORMALIZED DEL TMP X tCp 


DELTMP--T(t) - T(T=0) (Of ) 


CPj 108 


FIGURE 6-18 

REPRODUCIBILITY DATA 
6IN.TANK-STRAT TEMP PROFILES (r = 1.0) 
TESTS 8G NO. lOS AND NO. 16S (L -t- U) HEATING 
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TANK HEIGHT 


CR 1081 



FIGURE 6-19 

REPRODUCIBILITY DATA ; 

6 IN. TANK -STRATTEMP PROFILES (r= 1.0)' 
TESTS 8G NO. 1, NO. 5, AND NO. 8S (L) HEATING 


LIQUID LEVEL 


NO. 8S (L) 87.5% FILL 
:i^;8G NO. 5(L) 87.5% FILL 
iHif;8G NO. 1 IL) 87.5% FILL 
a =600 BTU/HR Ft2 




DELTMP = T(t)-T|t=0) (°F) 


0.6 0.8 1,0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 

NORMALIZED DEL TMP X 10^ 
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NORMALIZED TANK HEIGHT 


CR 1081 


FIGURE 6-20 

RErRODUCIBILITY DATA 
6 IN. TANK - STRAT TEMP PROFILES (r = 1.0) 
TESTS 8G NO. 4 AND NO. 6 (L + U) HEATING 
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LIZEDTANK HEIGHT 


CR 108-1 


FIGURE 6-21 

REPRODUCIBILITY DATA 
12 IN. TANK - STRAT TEMP PROFILES {T = 1.0) 
TESTS 1G NO. 30 AND NO. 25 (L) HEATING 



77 





CR 108-1 


FIGURE 6-22 

REPRODUCIBILITY DATA 
12 IN. TANK - STRAT TEMP PROFILES (r= 1.0) 
TESTS 1G NO. 31 AND NO. 17 |L + U) HEATING 



6 

12 


18 24 30 ' 42 

DELTMP = T(t) -T(T = 0) (°F) 

46 

60 


5 6 

8 10 


14 16 18 20 22 24^ 33 34 

NORMALIZED DEL TMP X 10^ 

36 

38 40 
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NORMALIZED TANK HEIGHT 


CR 108-1 


FIGURE 6-24 

REPRODUCIBILITY DATA 
12 IN. TANK -STRAT TEMP PROFILES (r= 1.0) 
TESTS 1G NO. 34 AND NO. 35 (L + U) HEATING 



0 4 10 


20 30 4 0 50 

DELTMP = T It) - T I t = 0) (®F) 


mmmmammmmmmmmmmmmmmm 

0.6 1.0 2.0 3.0 4.0 5.0 

NORMALIZED DEL TMP X 103 


6.0 


7.0 
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CR 108-1 


FIGURE 6-25 

REPRODUCIBILITY DATA 
12 IN. TANK -STRAT TEMP PROFILES (r = 1.0) 
TESTS 1G NO. 20 AND NO. 22 (L) HEATING 



6 8 10 12 14 16 

DEL TMP = T (t) - T (T = 0) (Op) 

0.8 1.0 1.2 1.4 1.6 1.8 2.0 

NORMALIZED DEL TMP X 103 






FIGURE 6-26 

REPRODUCIBILITY DATA 
12 IN. TANK -STRAT TEMP PROFILES (r= 1.0) 
TESTS 1G NO. 23 AND NO. 28 (L + U) HEATING 
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NORMALIZE 


CR 108-1 



FIGURE 6-27 

REPRODUCIBILITY DATA 
12 IN. TANK - STRAT TEMP PROFILES (r= 1.0) 
TESTS 1G NO. 26 AND NO. 15 (L) HEATING 
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NORMALIZED TANK HEIGHT 



FIGURE 6-28 

REPRODUCIBILITY DATA 
18 IN. TANK - STRAT TEMP PROFILES (r = 1.0) 
TESTS 1G NO. BS AND NO. D (L) HEATING 
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THERMOCOUPLE NO. 


CR 108-1 


FIGURE 6-29 

REPRODUCIBILITY DATA 
18 IN. TANK - STRAT TEMP PROFILES 
TESTS 1G NO. (F-6-18), NO. BS, AND NO. D 







50% LIQ LEVEL 
Gr* » lo"'® 

BTU 

q " LIQ HTG ONLY - 745—' 

^ HR FT2 

t >9 MINS 

— 1G - NO. F-6-18 (2-16-71) 

1G NO. D (7-18-72) 

1G f'JO. BS (7-18-72) 













LIQUID LEVEL 





0.25 0.5 

1 

2 3 ^ 

1 5 



NQRMALIZED4TX 10^ 



6 8 
ilT-TT-T.r.0‘°'=’ 
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FIGURE 6-30 cr los i 

REPRODUCIBILITY DATA 
181N. TANK-STRAT TEMP PROFILES (T= 1.0) 



Ar^rT-rr.oi°f) 
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THERMOCOUPLE NO. 


CR 108-1 


37 

38 

39 

40 
41-47 

48 

49 

50 
51-57 

58 

59 

60 


FIGURE 6-31 

REPRODUCIBILITY DATA 
12 IN. TANK -STRAT TEMP PROFILES (T= 1.0) 
TESTS 1G NO. A, NO. 30, AND NO. 25 






BULK LIQUID TEMPERATURE REPRODUCIBILITY CURVES 


FIGURES 6-32 THROUGH 6-45 


NOMENCLATURE 


q'^ Measured heat heat-flux input, Btu/hr-ft^ 

"s" Symbol in test number indicates no 

destratification heating 

(L) Liquid heating only 

(L+U) Liquid and ullage heating 

Liquid Freon 113 (PCA) 

T = (t/6) where x - Normalized time 

t = Elapsed test time 

6= Stratification test time 

6-In. Tank. . . . Title refers to tank diameter 
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CR 108-1 


FIGURE 6-32 

REPRODUCIBILITY DATA 
6 IN. TANK - BULK LIQUID TEMP HISTORIES 
TESTS 8G NO. 2 AND NO. 12S 


00 

(O 



0.8 1.0 1.2 
NORMALIZED TIME 









CR 108-1 


FIGURE 6-34 

REPRODUCIBILITY DATA 
6 IN. TANK - BULK LIQUID TEMP HISTORIES 
TESTS 8G NO. 9S AND NO. IBS 


(O 



0.8 1.0 1.2 
NORMALIZED TIME 







700BTU/HR FT 


STRATIFICATION 


DESTRAliFICATION 


FIGURE 6-35 

REPRODUCIBILITY DATA 
6 IN. TANK - BULK LIQUID TEMP HISTORIES 
TESTS 8G NO. 10S AND NO. 16S 


CR 108-i 


0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 

NORMALIZED TIME 
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FIGURE 6-37 

REPRODUCIBILITY DATA 
6 IN. TANK - BULK LIQUID TEMP HISTORIES 
TESTS 8G NO. 4 AND NO. 6 



^*H5gHtgr:3H 

8G NO. 4(L + U)87.5% FILL tig pn i jiff j? 
8G NO. 6 {L + U) 87.5% FILL JigiiiHpl g 

600 BTU/HR ^ ^ i iifeli i 








iliilli lii IMliiliiiii 


•::*nr:nrn:i:T:iH:^liHnn^fc^riHrhrrr5i^H~iHr!&rrH 


iiip 


f!:-; 

iiiii 













iui 


DESTRATIFICATION 






NORMALIZED TIME 






t;;:: 


. T . . - . - ■ I ^ - - » » T I I 

. T ■ ■ - i j j j I ..jixjjjjl 







CO 

CJl 



FIGURE 6-38 

REPRODUCIBILITY DATA 
12 IN. TANK - BULK LIQUID TEMP HISTORIES 
TESTS 1G NO. 30 AND NO. 25 


1G NO. 30 (L) 87.5% FILL 
1G NO. 25 (L) 87.5% FILL 
= 500BTU/HR Ft2 


STRATIFICATION 


destratification 


NORMALIZED TIME 



TEMPERATURE TF) 



tiiilS 


1G NO, 31 (L+U) 87.6% FILL 
16 NO. 17 (L+U) 87.5% FILL 


500 BTU/HR FT 


STRATIFICATION 


DESTRATIFICATION 


FIGURE 6-39 

REPRODUCIBILITY DATA 
12 IN. TANK - BULK LIQUID TEMP HISTORIES 
TESTS 1G NO. 31 AND NO. 17 


CR 106-I 


NORMALIZED TIME 







•1GN0. 33 (U 50% FILL 
"IG NO. 32 (L) 50%FILL 
350BTU/HR FT2 


STRATIFICATION 


DESTRATIFI CATION 


FIGURE 6-40 

REPRODUCIBILITY DATA 
12 IN. TANK - BULK LIQUID TEMP HISTORIES 
TESTS IG NO. 33 AND NO. 32 


CR 108-1 


0.8 1.0 1.2 

NORMALIZED TIME 




TEMPERATURE ( 
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STRATIFICATION 


DESTRATIFICATION 
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FIGURE 6-42 

REPRODUCIBILITY DATA 
12 - IN . TANK - BULK LIQUID TEMP HISTORIES 
TESTS 1G NO . 20 AND NO . 22 





-•1G NO. i) (L) 87.5% FILL |t 
-*1G NO. 22 (L) 87.5% FILL if 

= 300BTU/HR Ft 2 in 
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FIGURE 6-43 

REPRODUCIBILITY DATA 
12 IN. TANK - BULK LIQUID TEMP HISTORIES 
TESTS IG NO. 23 AND NO. 28 





OGNO. 23 IL + U) 87.5% FILL 
■•IG NO. 28IL + U) 87.5% FILL 
300 BTU/HR F t2 II 




^1 


lill 

;ur 






I li li ‘ ' 11 AT I idl iiiii 

lip 





-j:: H 














iilBiliilliiii 

I Dgl|VTf FjC^IC^ 


tH! fei i:U'^ 




tin :::; 


ilil 




iiil 

Hi: Hi: hi: 
Hr: ♦::: rrrrttr 

:t:t :ii: 

’“•i: :::: THi 

Eilj -p j--| 


NORMALIZED TIME 



FIGURE 6-44 

REPRODUCIBILITY DATA 
12 IN. TANK - BULK LIQUID TEMP HISTORIES 
TESTS 1G NO. 26 AND NO. IS 


CR 108-1 
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FIGURE 6-45 

REPRODUCIBILITY DATA 
18 IN. TANK - BULK LIQUID TEMP HISTORIES 

TESTS IG NO. BS AND NO. D 


-• 1GNO. BS(L)50%FILL 
-• IG NO. DID 50% FILL 
748 BTU/HR Fj2 


NO DESTRAT HTG 


STRATIFICATION 


NORMALIZED TIME 
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Table 6-2 


LIQUID BULK TEMPERATURE DEVIATION (PERCENT) 
(REPRODUCIBILITY) 


6- 

. 

Inch-Diameter Tank 

Grg 

Tb Deviation 
(Percent) 


50% Fill 

Liquid Heating 

10l4 

0.68 



Liquid- Ullage Heating 

10l4 

0.09 


87.5% Fill 

Liquid Heating. 

7 X 1014 

0.22 




1015 

0.13 



Liquid-Ullage Heating 

7 X 1014 

0.69 




1015 

0.21 

12 

-Inch- Diameter Tank 




50% Fill 

Liquid Heating 

10 14 

0.17 



Liquid-Ullage Heating 

10 14 

0.43 


87.5% Fill 

Liquid Heating 

10 14 

0.3 




7 X 10 14 

0.35 




10l5 

0.04 



Liquid- Ullage Heating - 

7 X 10l4 

1.03 




1015 

0.18 

18 

-Inch- Diameter Tank 




50% Fill 

Liquid Heating 

1015 

0.21 


remain apart over the entire test period. However, in all oases the curves 
remain nearly parallfil emphasizing that the rate of increase of the bulk tern 
perature for each test is nearly identical. 


The temperature daviatioh values are shown to be less than about; 1 percent as 
indicated in Table 6-2, The liquid- bulk temperature represents the integrated 
value of the respective temperature profile curve and shows that for each 
repeated test an equal amount of thermal energy has been transfered to the 





t03 









































liquid. It should be noted, however, that for some comparisons, the test 
conditions for t >1.0 were no longer identical and comparisons made beyond 
this time have ho meaning. This was true since in some tests the heaters 
were not turned off during the destratification period. 

* 

Based upon the reproducibility shown in Jrigure 6-15 through 6-45 and the 
deviation values shown in Table 6-Z, it is concluded that ahigh degree of reli- 
ability in reproducing the temperature values has been achieved. Thus, these 
data have demonstrated that each temperature profile ’'signature, ” or 
temperature anomaly, at each respective position has been reliably 
reproduced for each set of test conditions for each tank. This particular 
signature trait reproducibility is ah important factor in substantiating the 
data reliability where in a sense every wrinkle in the temperature profile 
can be repeated for a given tank. 


6. 3 EXPERIMENTAL MODIFIED GRASHOF NUMBER DEFINITIONS 


6. 3. 1 Wall Heat Flux Variation Effects on Modified Grashof 
Number Definitions 

The Modified Grashof number was defined in Section 2 as gpq''L'^/kv^ where 
q“ was listed as the wall heat flux. However, during the course of the 
experimental program, it became necessary to delineate the difference 
between the heat flux which is generated within the heater blanket, q^, and 
the heat flux to the liquid, qT. Hence, two different Modified Grashof num- 
berS, Gr, 


4 > 2 

L /kv ), are used in discussing 


•jj ^ ) 3-Ud Gr^ 

the data. The experimental value of Gr^for any test is constant with time 
since qj^, is kept constant. However Gr* varies with time since q^b will 
differ from q^ by the rate of energy storage in the tank walls and the heater 
blanket. In fact, both Gr^ and the Interface number {h^gk^/q^V Dc^ ) change 
by more than an order of magnitude between t = 0 and t = 1. 0. ^ 


The term Gr^ is the scaling term which must be used when treating the total 
tank s y s t em. , Her e, the tank wall s , the in sulatio n, and the he ate r s , in addition 
to the liquid and ullage, make up the total system. Since the three tanks are 
s caled thermo dynamically and geometrically and have identical; material, it is 



only necessary that "which reflects the boundary heating condition be the 
same for each tank. When this dimensionless group is kept the same for two 
separate tanks along with the Fourier and Interface numbers, then the two 
tanks are considered scaled to one another. 

It is important to realize, however, that if the liquid and ullage are treated as 
the thermod-ynamic system, that Gr^ and not Gr^ will have meaning with 
respect to this system. This can be explained by examining the equation 

= q'^ + q'g (6-3) 

where qh equals the total heat-flux output from the electrical heater 
blankets, q^' equals the heat-flux transfered to the liquid at any time, t, and 
q’* equals the instantaneous rate of accumulation or storage of thermal 
energy per unit area of tank wall, insulation, and heaters. Since q^' depends 
on the tank walls having an inside surface temperature greater than that of the 
contained fluid, the major portion of q^ in the early part of the test period 
will be accounted for by q'J. 

Figures 6-46 and 6-47 are typical examples of the trend of the measured 

Sjt ■ ■ ' 

values of Gr'p' with respect to Grjj as a function of time for this study. As can 
be observed, the ratios of GrJ/Gr^ (which is also equal to q^'/q^) are less 
than one over the test period although their values are increasing with time. 

It is likely that this ratio will eventually reach a constant value, but the data 
are insufficient to show this value. Figures 6-48 through 6-Bl show typical 
curves of Gr-g and Gr^ . It should be noted that the value of Grjj in these 
figures increases slightly with time even though is kept constant; This 
change reflects the liquid property changes-with temperature that occur 
during the course of a test. However, the initial value of Grjj (i. e. , at t = 0,0) 
is the one referred to in this report. The Gr^ and Gr^ values at 1,0 for 
each test are shown in Appendix A. Volume II of this final report includes the 
Grjj and GrJ curves as a function of the normalized tirne: for each experiment. 

In summary, there are three main points that have been made. These are: 

A. The independent scaling parameter Gr.g is' the proper Modified 
Grashof number for scaling a hank system. 



FRACTION OF HEAT TRANSFERRED TO LIQUID 
STRATIFICATION TESTS 


FIGURES 6-46 AND 6-47 


NOMENC LATURE 


Gr”^ . ... Modified Grasliof number 

q” * • ♦ • Measured heater heat flux at source 

'H 

q'' Determined liquid heat fluXj 

experimentally absorbed 

(L-tU) . . . Liquid and ullage heating 

(L) Liquid heating 
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FIGURE 6-46 

FRACTION OF HEAT TRANSFERRED TO LIQUID-STRATIFICATION TESTS 


* X 
u 
CS 

s 

II 

cr 



0.4 0.6 

NORMALIZED TIME 


FIGURE 6-47 

FRACTION OF HEAT TRANSFERRED TO LIQUID-STRATIFICATION TESTS 


CR 108-1 



0.4 0.6 

NORMALIZED TIME 
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LIQUID MODIFIED -GRASHOF NUMBER HISTORY CURVES 
STRATIFICATION TESTS 

FIGURES 6-48 THROUGH 6-51 

NOMENC LATURE 

(3. . . Coefficient of tliermal expansion of fluid 

p. . . Fluid density 

[X. . . Fluid viscosity 

k. . , Fluid thermal conductivity 

T. . . Fluid temperature 
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FRECN 113 {PCAl MODIFIED GRASHOF NO. (Gr^l 
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FIGURE 6-50 

12 IN. TANK 1G NO. 34 LIQUID MODIFIED GRASHOF NO. HISTOmES-STRAT 


= 350 BTU/HR FT^ JL + U - HTGi 


q = asu 
50% FILL 


1 


K-{Gri) 


0 0^ 

' ■ « ^ 
0.4 0.6 0,8 1.0 


normalized TIME 


1 -THEORETICALLY BASED ON MEASURED q” WATTAGE INPUT 

H 

X- EXPERIMENTALLY BASED ON DETERMINED q" ; ' 

L= fluid depth 
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FREON 113 (PCAl MODIFIED GRASHOFNO. (Gr'J 


CR 108-1 

FIGURE 6-51 

■{2 !N. TANK1G N0;33 LIQUID MODIFIED GRASHOF NO. HISTORIES-STRAT 



■ ■ ■. /. . 

= f \ 1- theoretical based ON MEASURED WATTAGE INPUT 

V (c ’ ] X- EXPERIMENTAL BASED ON MEASURED q" ' 

L= fluid DEPTH 





B. The Modified Grashof number, is a dependent parameter 

obtained experimentally, that might be useful wheu fluid and ullage 
are considered as a system. 

C. The value of Gr^ 'will range from zero to a value less than 
during the stratification period. 

6. 4 CRITERIA FOR GAGING DEGREE OF SCALING 

The following subsections compare the dimensionless pressure and dimension- 
less temperature data for scaled tests. For such scaled tests, the Grjg;, Fo, 
and I are equal for the two tanks. Hence, for percent scaling it would be 
expected that the dimensionless pressure and temperature values be the same 
for these scaled jtests. However, no established criteria to measure the deg- 
gree of success in scaling are available. Thus, the criteria used here were 
selected to represent both the integrated or volumetric heat transfer effect 
as well as a pointwise effect. The requirement that the bulk temperatures 

= TbeM (5-4) 

is at best a necessary requirement for scaling and reflects the integrated heat 
transfer effect. On the other hand, the temperature profile requirement at a 
point 

TA (x,y;i,T) = Tb (x, y,z, t) (6-5) 

/ ■ 

if met everywhere in both liquid and the ullage would likely provide a necessary 
and sufficient condition for complete scaling. The pressure requirement 

Pa = Pb ,(6-6) 

I • ' ' 

may lie somewhere in between the above two criteria because in effect it is an 
integrated value and yet is coupled to the liquid surface temperature. The 
extent to which these three conditions are met wi ’1 Gharacterize the degree 
of scaling that was accomplished in two scaled tank tests, 

A quantitative guide for the degree, of scaling can probably best he considered 
in terms of the percent deviation between the values of p, Tb, and T in a set 
of scaled tests. The percent deviation is presented for each term, p, Tb, and 
T, in Subsections 6.5 and 6. 6. 



6. 5 PRESSURE HISTORY SCALING RESULTS 
6. S, 1 latyoductipn 

Tile pres.sui'e history scaling curves are platted in both the dimensional 

j 

pressure^ p, and the dimensionless pressure, p, on the same graph. This 
can be done since p and p only differ by the constant 6/D) which is 
divided into p to obtain p. This term is constant for any two scaled tests 
using the same liquid. Pressure measurements were obtained both with a 
pressure transducer and a pressure gage. Subsequent calculations used the 
arithmatic average of these two experimental pressure values. However, 
only the pressure transducer was used in the high-g tests. 

The stratification portion of each experiment corresponds to values of r less 
than or equal to 1.0. All scaled tests considered here, with liquid heating 
only, were scaled reasonably well in this regime. The liquid- ullage heating 
tests did not compare as favorably, 

6.5.2 Scaling Evaluation 

Figures 6-52 through 6-61 show the pressure as a function of time for sets 
of scaled tests covering the range -of tank configurations discussed earlier. 

The most general observation that can be made is that the matched sets of 
curves agree well for liquid- only heating while the agreement is not as good 
for the liquid and ullage heating cases. The extent to which scaling was 
achieved may be judged on the basis that the dimensionless pressure curves 
definitely indicate the same general trends and values. This should be judged 
taking into account that the two tanks may have different operating conditions. 
As an example, when the 6-inch tank experiences 27 g*s with 2, 230 Btu/hr-ft , 
the scaled equivalent is an 18-inch tank at 1 g with 748 Btu/hr-ft^. Thus, 
when the dimensionless pressure curves agree reasonable well with entirely 
differ ent .operating conditions, it is likely a r esult of scaling, 

One method of comparison uses the pressure difference values between t = 0 
and T - 1.0. Here, the rise in pressure for one tardc is compared With that of 



SCALING PRESSURE DATA 
TANK -ULLAGE -PRESSURE HISTORIES 


FIGURES 6-52 THROUGH 6-61 


NOMENCLATURE 

q” ...... . Measured heater lieat-£lux: source, Btu/hr-ft^ 

H 

Q ........ Stratification, test period 

D ...... . . Tank diameter 

Gr^. ...... Modified Grashpf number based upon q”H 

CHPRESS . . . Characteristic pressure, (q"jjO/D) 

p , . . . .... Tank- ullage measured pressure 

. Normalized pressure, (p/CHPRESS) 

6 in. -8G#2S . Test title (e. g. , 6-inch dia tank - 8 g's test 
, number 2; S indicates no heating during 
destratification testing) 



TANK-ULLAGE PRESSURE (PSIA) 
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FIGURE 6-52 

SCALING PRESSURE DATA - TANK PRESSURE HISTORY, 
6 IN. 8G2- 12 IN. 1G30 
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TANK-ULLAGE PRESSURE (PSIA) 


FIGURE 6-53 

SCALING PRESSURE DATA - TANK PRESSURE HISTORY 
6 IN. 8G9S- 12 IN. 1G33 


CR loe-i 



117 



















FIGURE 6-54 

SCALING PRESSURE DATA - TANK PRESSURE HISTORY 
6 IN. 8G1 - 12 IN. 1G20 


CR 108-1 














TANK-ULLAGE PRESSURE (PSIA) 


FIGURE 6-55 

SCALING PRESSURE DATA - TANK PRESSURE HISTORY, 
6 IN.27G1S- 18IN. IGBS 


CR 108-1 



NORMALIZED TIME 


119 







TANK-ULLAGE PRESSURE (PSIA) 


CR 106-1 


FIGURE 6-56 

SCALING PRESSURE DATA - TANK PRESSURE HISTORY, ' 
6IN.27G3S-18IN. 1GIS 
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TANK-ULLAGE PRESSURE (PSIA) 


FIGURE 6-57 

SCALING PRESSURE DATA -TANK PRESSURE HISTORY, 
6 IN. 8G14 - 12 IN. 1G26 AND 15 
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TANK-ULLAGE PRESSURE (PSIA) 
















TANK-ULLAGE PRESSURE {PSIAI 


CRIOe-l 


FIGURE 6-59 

SCALING PRESSURE DATA -TANK PRESSURE HISTORY, 
6 IN. 8G10S AND 16S - 12 IN. 1G34 
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tlxe other scaled tank and a percent difference is obtained by dividing this 
difference by the average initial pressure of both tanks, i.e., 

^Pl~Z 

p. 

lavg 

This method was also used to determine the deviation in reproducibility as 
shown in Table 6-1. Table 6-3 summarizes the percent deviation values for 
the scaled tests. Here again the best agreement was with the liquid-only 
heating tests as indicated in the qualitative evaluation where the percent 
deviation was 5.1 percent or less. This would suggest in a quantitative 
manner that liquid- only heating tests were reasonably well scaled with 
respect to pressure. 

The liquid-ullage heating tests did not scale as well. It is to be noted that 
the error band when testing for repeatability was up to about 8. 4 percent 
between identical tests with the same tank. The measured difference for the 
: scaled tests were between l6 to about 28 percent for the 87. 5 percent fill 
cases. Thus, even if a likely spread in data of 10 percent is accounted for, 
the additional 6 to 18 percent difference is likely due to differences in the 
heat transfer phenomena not accounted for in. the scaling parameters. 

6. 6 BULK TEMPERATURES” TEMPERATURE PROFILES-SC ALING 
RESULTS 

6. 6. 1 Intro due tip h 

Temperature data are presented in bulk- liquid temperature-time curves, and 
temperature profiles at normalized time t = 1.0 for each test in both the 
dimensional and nondimertsional forms. The temperature is nondimensionalized 
by dividing it by the characteristic temperature (q^ D/k^) which is constant 
for any set of scaled tests. Temperature profile curves at several stratifica- 
tion test times as well as and the liquid, ullage, and tank structure tempera- 
ture data for each test are presented in Volume II. 

The bulk liquid temperaturo histories for the stratification scaling tests 
scaled to within 3 percent. However, a greater deviation was found when 


- Pii 


1/2 (Pi. + P2j) 


P2f " PZj. 
1/2 (Pi. ¥ P2.)‘ 


(6-7) 



Table 6-3 

PRESSURE DEVIATION (PERCENT 
(SCALING TESTS) 


8G Tests — 6-In.- 
IG Tests - 12-In 

Dia Tank/ 
.-Dia Tank 

G?h 

Deviation 

(Percent) 


50% Fill 

Liquid Heating 

10l4 

4.4 



Liquid-Ullage Heating 

10 14 

16,2 


87.5% Fill 

Liquid Heating 

10l4 

5.1 




7 X 10l4 

1.8 




1015 

2.0 



Liquid- Ullage Heating 

10l4 

16.0 




7 X 10l4 

27.8 



• 

iqIS 

23.0 

27G Tests — 6-In,Dia Tank/ 
IG Tests ™18-In.-Dia Tank 




50% Fill 

.... 

Liquid Heating 

10l5 

4.9 


87.5% Fill 

Liquid Heating 

ioi6 

5.2 


making a pbintwise comparison of the temperature profiles for the scaled 
tests. ■ . ‘ ■ 

6.6,2 Temperature Scaling Results 

6. 6. 2.1 Bulk Temperatures 

Nondimensional aud dimensional bulk liquid (Freon 113) temperature history 
curves are shown for the stratification scaling tests in Figures 6-62 through 
6-71. The normalized bulk liquid temperature history curves indicate that the 
scaled tests agree 'where the curves match well for both the liquid and liquid- 
ullage heating conditions. The bulk temperature reflects the integrated tem- 
perature effect of the heat absorbed by the liquid and is a measure of the 
taiik-system scaling concept studied in this program. Hence, qualitatively, 
the curves indicate that the total energy transferred to the liquid has been 
scaled reasonably well. 
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SCALING TEST DATA 
BULK LIQUID TEMPERATURE HISTORIES 
(FREON 113 (PCA)) 


FIGURES 6-62 THROUGH 6-71 


NOMENC LATU RE 


Gtjj Modified Grashof number based on input heater 

heat-flux sov VC e, (g L'^/Tc^v2) 

g. .•••••• • Gravitational constant 

• • Liquid coefficient of thermal expansion, 

Liquid density 

T Measured temperature (°F) 

q" . Measured input heat heat flux source, Btu/hr-ft^ 

N 

L Liquid depth 

kjg Liquid thermal conductivity 

Liquid kinematic viscosity, {H-i/p^) 

|j,^ ........ Liquid dynamic viscosity 

D, Tank diameter 

CHTMP .... Characteristic temperature, (q|^ 

T. Normalized temperature, (T/CHTMP) 

T. . . Normalized time, (t/G) 

t ........ . Elapsed test time ^ 

0 Stratification test time 

6in.-8Grrl2S . Test title (e. g. , 6- inch diameter tank- 8 g's 

test number 12; "S” indicates no heating during 
de stratification testing) 


) 



TEMPERATURE (°F) 
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FIGURE 6 62 

SCALING tests DATA - BULK LIQUID TEMP HISTORIES 
6 IN. 8G2 AND 12S - 12 IN. 1G30 AND 25 


LIQUID HEATING 
87.5% FILL 
Gru=10l5 


l^OPfSTRAT^H^Ti^Q 
I.-8GNO. t2S r T 


STRATIFICATION 


DESTRATIFICATION 


normalized TIME 
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FIGURE 6-63 

SCALING TESTS DATA - BULK LIQUID TEMP HISTORIES, 
6 IN. 8G3 AND IIS - 12 IN. 1G31 AND 17 


LIQUID + ULLAGE HEATING 
87.5% FILL 


NO DESTRAT HEATING 


STRATIFICATION 


DESTRATIFiCATION 


NORMALIZED TIME 









FIGURE 6-64 

SCALING TESTS DATA - BULK LIQUID TEMP HISTORIES, 
6 IN. 8G9S AND 15S - 12 IN. 1G33 AND 32 
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50%FILL .. 
Gr^io’^ 


STRATIFICATION 
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NORMALIZED TIME 
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FIGURE 6-65 

SCALING TESTS DATA - BULK LIQUID TEMP HISTORIES, 
6 IN.8G10S AND 16S- 12 IN. 1G34AND35 


LIQUID + ULLAGE HEATING 
50% FILL 
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DESTRATIFICATION 
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FIGURE 6-66 

SCALING TESTS DATA - BULK LIQUID TEMP HISTORIES, 

6IN.8G1,5AND8S- 12 IN. 1G22AND20 
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LIQUID + ULLAGE HEATING 
87.5% FILL . 

Gri=7X,0'" iteiMil 
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FIGURE 6-67 

SCALING TESTS DATA - BULK LIQUID TEMP HISTORIES, 
6 IN. 8G4 AND 6-12 IN. 1G23 AND 28 
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FIGURE 6-68 

SCALING TESTS DATA - BULK LIQUID TEMP HISTORIES, 
6 IN. 27G1S - 18 IN. 1GBS AND D 
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FIGURE 6-70 

SCALING TESTS DATA - BULK LIQUID TEMP HISTORIES, 
6 IN. 8G14- 12 IN. 1G26 AND 15 


LIQUID HEATING 


NORMALIZED TIME 
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FIGURE 6-71 

SCALING TESTS DATA - BULK LIQUID TEMP HISTORIES 
6 IN.8G13- 12 IN. 1G27 


NORMALIZED TIME 


A quantitative look at the degree of scaling or matching of the bulk 
temperature data is shown in Table 6-4. This table shows that the bulk liquid 
temperature curves scale to within a 2. 5 percent deviation. The basis for 
these estimates is discussed in Subsection 6. 2. 2 (Equation 6-1), The differ- 
ences expressed in these values and in Figures 6-62 through 6-71 may be 
directly related to the measured Modified Grashof number variations 
discussed in Subsection 6. 3. 1. 


6.6. 2. 2 Temperature Profiles 

Although the bulk-liquid temperature data show that the scaling criterion 
Tt)^(T) = Tbg (t) is met to within a 2. 5 percent deviation, th^ additional cri- 
terion of T^ (S, y, z, t) = Tg (5c, y, z, t) must also satisfied. The delta tempera- 
ture and normalized delta temperature profiles for the scaled tests at the end 
of the stratification period (t = i. 0) are compared below. 


Table 6-4 

MAXIMUM BULK TEMPERATURE DEVIATIONS (PERCENT) 

(SCALING TESTS) 


8G Tests“6-In.-Dia Tank/ 
IG Tests — I2~Tn.Dia Tank 

, 

Deviation 

(Percent) 


50% Fill* 

Liquid Heating 

1014 

0.8 



Liquid- Ullage Heating 

10^^ 

0.5 


87.5% Fill 

. Liquid Heating 

ioi4 

1,4 . 




7 X IQI^ 

1.1 




10l5 

0.9 



Liquid- Ullage Heating’ 

I0l4 

1.4 




7 X 10l4 

2.5; 




1015 

0.9 

^ 27G Tests — 6-In.-Dia Tank/ 
IG Tests — 18-In.-Dia Tank 




50% Fill 

Liquid Heating 

1015 

■' L? .. 


87.5% Fill 

Liquid Heating 

10l5 

0.9 
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These normalized delta temperature proliles at t = 1,0 for scaled tests are 
shown in Figures 6-72 through 6-81. The data show that there are differences 
in the normalized profiles for the scaled tests. Because of the high degree of 
repeatability of the liquid- bulk temperature as well as the point-by- point 
profile: agreement^ it is likely that these differences are largely due to real 
thermodynamic and fluid dynamic differences in the separate tanks. It should 
be noted, however, that the difference in temperature profiles of a scaled 
pair of tank tests was found to be greater than that between profiles for the 
same tank even when Gr^ differed by about an order of magnitude. This is 
illustrated in Figures 6-72 through 6-81 which show that the axial temperature 
gradient for the 6-inch tank at the 8-g and 27-g levels exceeds the 12-inch and 
18- inch tank gradients at 1 g for each scaled set of tests. This was found to be 
true for both 50 and 87.5 percent fill levels. Furthermore, Figure 6-82 shows 
that for the same 6-inch tank and for values of Grj^ from 10^"^ to 7 x that 

the temperature gradient difference was significantly less than that between 
scaled tests. This would seem to imply that there is an inherent difference in 
each of the tanks which may or may not be related to the scaling effect. 
Unfortunately insufficient data are available to pinpoint the cause of the tank- 
to -tank variations although the fact that the heater blankets were all of the 
same thickness may be a contributing factor. 


A quantitative assessment of the variations in the liquid- scaling temperature 
profiles was made by determining the percent deviation of the liquid- 
temperature profile slopes for the scaled tests. This was done by first 
determining the area under the profile curves. Then an average stratification 
temperature profile slope was determined for each curve. The percent 
deviation for each set of scaled data was then determined by dividing the 
difference in the liquid profile slopes by the mean slope for the scaled tests 
as shown: 


PROFILE DEVIATION =■ 


. 

CiTN . 

/0T\ 

1 

2 

(f) 

+ . 

j'm 


(6-8) 


Hence, Table 6-5 shows the profile percent deviations for the scaled 
experiments. These values are a relative measure of the deviation of the 
point-by-point liquid-temperature scaling requirement. 
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SCALING TEST DATA 

STRATIFICATION TEMPERATURE PROFILES (t = 1.0) 


FIGURES 6-7Z THROUGH 6-81 


NOMENCLATURE 


g4, 


g- 

T. 

q” 

L. 


D. 
CHTM 
AT 
AT 

T . 

t . 

0 . 


Modified Grashof number based on input heater 
heat flux source, {gp^ q^ L^/k^v^2) 

Gravitational constant ^ 

1 

Liquid coefficient of thermal expansion, (*^ ) 

Liquid density 

Measured temperature (°F) 

Measured input heater heat flux source, Btu/hr-ft^ 
Liquid depth 

Liquid thermal conductivity 
Liquid kinematic viscosity, 

Liquid dynamic viscosity 
Tank diameter 

Characteristic temperature, {q, 

X-t 

Del Tmp, T (t) - T (x = 0) 

Normalized Del Tmp, (AT/CHTMP) 

Normalized time, (t/0) 

Elapsed test tune 
Stratification test time 


D/kp 


6 in.-8G#l2S. Test title (e. g. , 6- inch diameter - 8 g^s 
test number 12; *'S’* indicates no heating 
during de stratification testing) 
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1^87.6% LIQUID LEVEL 
itlFLIQUID + ULLAGE HEATING 
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50%LIQUID LEVEL 
LIQUID HEATING 
* 14 

Gr u = 10 
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NORMALIZED TANK HEIGHT 
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NORMALIZED TANK HEIGHT 


CR 108-1 



87.5% LIQUID LEVEL 
LIQUID + ULLAGE HEATING 

Gru - 7 X 


r^lGURE 6-77 

SCALING TESTS DATA - STRATIFICATION TEMP PROFILES (r = 1.0) 
6 IN. 8G4 AND 6 - 12 IN. 1G23 AND 28 


1! 3 4 

NORMALIZED DEL TMP X 10^ 


20 

DELTMP = T(r)-T|r=0) (°F| 
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NORMALIZED TANK HEIGHT 
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NORMALIZED TANK HEIGHT 
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NORMALIZED TANK HEIGHT 



87.5% LIQUID LEVEL 
LIQUID + ULLAGE HEATING 
. 14 

Gr u ■ 10 
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6-INCH DIAMETER TANK 
STRAT TEMPERATURE PROFILES 

(t= i.o)-"G" dependence 


FIGURE 6-82 


NOMENC LATURE 


Measured input heater heat-flux source, 

“ Btu/hr-ft2 

T Measured temperature (°F) 

AT Del Tmp, T (t) - T (t = 0) 

'AT Normalized Del Tmp, (AT/CHTMP) 

CHTMP ... Characteristic temperature, (q^ 

D. ....... . Tank diameter 

Liquid thermal conductivity 
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NORMALIZED TANK HEIGHT 



i 87.5% LIQUID LEVEL 
LIQUID HEATING 
! qi!!/ =600BTU/HR FT 


DELTMP 


0.8 1 


0.7H 


0.6H 


C.bH 
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Table 6-5 


LIQUID STRATIFICATION TEMPERATURE PROFILE 
DEVIATION (PERCENT) 

(SCALING TESTS) 


8G Tests — 6-In. - 
IG Tests-12-In. 

Dia Tank/ 
-Dia Tank 

Gr^ 

Deviation 

(Percent) 


50% Fill 

Liquid Heating 

10^4 

39.1 



Liquid-Ullage Heating 

10l4 

32.4 


87.5% Fill 

Liquid Heating 

10l4 

29.3 



■ 

7 X 10l4 

35.3 




10^5 

32.8 



Liquid-Ullage Heating 

10l4 

24.4 




7 X 10^4 

29.8 




10l5 

21.4 

Z7G Tests — 6-In. 

-Dia Tank/ 



IG Tests — 1 8-In.- Dia Tank 



■[ 

50% Fill 

Liquid Heating 

10l5 

18.4 


87.5% Fill 

Liquid Heating 

ioi6 

10.1 


Therefore, it can be concluded that the test system has been reasonably 
scaled to the bulk temperature deviation values which are approximately 
3 percent (Table 6-4); whereas the scaling of the stratified liquid temperature 
profile (gradient) results in profile deviation values from 10 to 39 percent 
(Table 6-5). 

However, it is to be noted that the liquid surface temperature differences, for 
the liquid heating only scaling tests, were within a range of about 3°F. Since 
the tank vapor pressure is indirectly coupled to the liquid surface temperature, 
it is interesting to note that the measured tank pressures for these tests were 
within about 1 psi of the saturated pressure corresponding to the measured 
Freon 113 liquid surface temperature. This was not true for those cases 
where the ullage was also heated. 


154 



Section 7 
CONCL.USIONS 

The experiments outlined in the matrix. Figure 5-1, have been conducted with 
three scaled cylind rical tanks, 6, 12, and 18 inches in diameter, with hemi- 
spherical domes having a tank length -to -diameter ratio of 3, These tests have 
been conducted with liquid (Freon 113) heating as well as with liquid and ullage 
heating, at both the 50 and 87. 5-percent filling capacities, with heat flux 
values from 45 to 2, 280 Btu/hr-ft^, Modified Grashof numbers (Gr^) from 
10^^ through 10^^, Fourier numbers from 10“^ through lO""^, and Interface 
numbers from 7. 54 x 10 “3 through 2, 42x 10"^. The experiments are com- 
pared in scaled sets by matching three dimensionless parameters (i. e. , Gr^^, 
Fo, and Interface number) obtained in Subsection 2. 5. The following con- 
clusions can be drawn from these experiments: 

A, The bulk-liquid temperature for the stratification scaling tests was 
scaled to within 2, 5 percent for all test configurations. The 
reproducibility data scatter of the bulk liquid temperature was less 
than 1 percent. 

B, The dimensionless tank pressure with liquid heating was scaled to 
within 5. 1 percent. During the reproducibility tests, data scatter 

for these cases varied between 0, 7 and 3. 5 percent with one exception. 

C, The dimensionless tank pressure of two scaled tanks with liquid 
and ullage heating was scaled in the range of 16 to 28 perce. . 

During the reproducibility tests, data scatter for this pressure and 
heating condition varied betvi'een 2 and 8.4 percent. 

D, The stratified liquid scaling profile deviation vaLues showed that 
scaling was achieved in the range of 10 to 39 percent, although the 
liquid surface temperature deviation values were about the same 
as the pressure deviation values for the liquid heating only cases. 

The profile reproducibilities were nearly identical to one another; 
hence, the deviations were shown to be neglieible. 
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In BUinmary, the data reproducibility for the stratification experiments was 
good and, as a result, provided an excellent basis for comparing scaled 
tests. It has also been shown that the scaling criteria for the bulk 
temperature in all experiments and the tank pressure with liquid heating only 
were satisfied with respect to the pi-groups, F^, and I, Thusj it is 

concluded that these dimensionless parameters can be used to scale when 
the quantities which reflect the integrated heat transfer effects, and p^j 
are of interest. In addition, the temperature deviation values for the liquid 
surface were small. However, in general, deviations in the pointwise tem- 
perature T(x, y, z, t) for the scaled tests have been found but the data are 
insufficient to determine whether or not these deviations are due to the 
inadequacy of the- pi groups or to a lack of scaling of the entire tank system 
with the constant heater thicknesses. 



Section 8 

RECOMMENDATIONS 

The following recommendations are made within the context that the prediction 
of tank pressures for heated propellant tank systems is necessary. This 
report indicated that certain dimensionless parameters were necessary to 
define a heated container of liquid when considering pressure and bulk 
temperature. Thus, certain steps will have to be taken if the present state 
of the art is to be extended towards the above goal of making pressure and 
temperature predictions. These steps will include the extension of the 
experimental range so that correlations can be made, comparing these data 
obtained here with existing experimental correlations, and a further examina- 
tion of the present data. 

Hence the recommended steps to be continued from this program are primarily 
in three phases; 

1. To evaluate and more fuUy understand the implications of the strati- 
fication (and destratification) data collected in this program. 

2. Further testing to expand the range of parameters both for determin- 
ing the important pi groups and correlating the tank pressure and 
stratified liquid temperature profiles with the important scaling 
parameters, 

3. Comparing data obtained here with existing experimental 
correlations. 

4. Using latest experimental correlations modify an existing or develop 
a new model to predict tank pressure and propellant temperatures. 
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Appendix A 

LIQUID AND ULLAGE PROPERTIES 


This appendix contains the test liquid (Freon 113 PCA) and ullage properties, 
and the characteristic temperature and pressures for each experimental test 
investigated. Also shown are the values for the three scaling parameters 
that were maintained constant for the scaling test in this study. These 
parameters are the Modified Grashof number (Gr^), the Fourier number 
(Fo), and the Interface number (I). This information is shown in Table A-1. 


Notes: 



Measured heater heat -flux source, Btu/hr-ft 


q 


II 

i 


Determined heat flux from heat absorbed by liquid, 
Btu/hr-ft^ 


T - 1 ......... Normalized time at end of stratification test. 

CHTMP Characteristic temperature, (q^ D/k ) 

CHPRESS .... Characteristic pressure, (q^ 9/D) 
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Table A- 1 


FREON 113 PCA LIQUID - ULLAGE PRC 


Teat Wo, 

6 in. ec It^s 

6 in. 8G HZ 

6 in. SC III 

6 in. ec i 

Characteristic Temperature (CHTMPJ (°F) 

8.S61 +03 

1, 367 + 04 

7. 595 + 03 

1. 140 + 03 

Characteriatic Pressure (CHPRESS) Ub^/ln^) 

4.311 +03 

3. 609 + 02 

3.96 + 03 

4,058 + 02 

LIQUU3: density (Ib/ft^) 

94.411 - 91. 656 

94. 3B2 - 93.681 

94. 387 - 93,539 

■ 

94, ^19 - 93. 1 
'i 

mass (lb) 

13.361 - 13. 359 

21.621 - 31.46 

21. 633 - 31.435 

31,61 - 31. 43 

speciftc heat (Btu/Ib 

0. 335 - 0. 337 

0, 335 - 0. 227 

0.235 - 0.337 

0, 325 - 0-227 

thermal conductivity {Btu./hr ft ^F] 

0. 0+ - 0. 019 

0.04 - 0,039 


0, 04 - D. 039 

1 viscoBlty {Ib/ft hr) 

1.302 - 1. 137 

J. 19fl - 1. 139 

1.198-1.127 

1. 194 - 1.13 

hinernatic viacosity (ft^/hr) 

0.013732 - 0.013139 

0.01369 - 0.01216 

0.012696 - 0.D12D49 

0,0U65 * 0,C 

Prandtl No. 

6. B34 - 6.6116 

6.819 - 6.6185 

6.83 - 6, 57B 

- 6. 5S*i 

■ 

Modified Graehof No. - { (t^^a ^ )a 

1. 222 + X4 

1, 343 + 15 

8. 317 + 14 

4 14 

j 

Modified Grushof No. = f {q”, t = 1, ), Gr^ 

6.413 + 13 

7. 8 + 14 

5. 33 + 14 

6+999 f ; 

. Interface No# , X 

0.03106 

0,03i8 

0.0363 

0.243 i 

ULLAGE; density (Ib/ft^) 

0.451 - 0. 531 

0.453 - 0, 546 

0,454 - 0. 559 

0.455 - 0. 56( 

mass (lb) 

0.0S9 - 0.07 

0.015 - 0.016 

0.015 - D.OIS 

0.015 - 0.01' 

specific licat (Btu/lb ’ F) 

0.158 - 0. 16 

0. 158 - 0. 16 

□. 158 - 0. 16 

0. 1S0 - 0,16! 

• i 

Test No. 

6 in. 27C fflS 

13 in. 1C (125 

13 in. IG 032 


Characteristic Temperature (CHTMP) (“F) 

3.663 +04 

1. 357 +04 

8.861 + 03 


Characteristic Pressure (CHPRESS) Ubj/ln.^) 

S. 195 +03 

3. 606 + 03 

4.094 + 03 

4. 053 +02 

■ J 

LIQUID: density (Ib/ft^) 

94. 139 - 93.009 

94. 33 - 93. 53 

93.42 -.93.39 

94,51 - 93. 3^ 
' 

mass (lb) 

21. 56S - 31.306 

173.85 - 171.41 

97.83 - 96.755 

173.63 - 171. 

spccinc heat (Btu/Ib *F] 

0. 226 - 0. 229 

0. 335 - 0. 227 

0.338 - 0.231 

0. 325 - 0. 33( 

thermal conductivity (Btu/hr ft *F) 

0.039 

0.040 - 0.039 

0,039 - 0.038 

0. 04 - 0, 039 

viscosity (!b/fb hr) 

1. 177 - 1.088 

T, 193.- 1. 137 

■ 19 - 1. 044 

1. 191 - i.ii: 

jj 

kinematic viscosity (ft^ /hr) 

0,0135 - 0.01169 

0.01364 - 0.01205 

U.0U97 - 0.0113 

; 

0.Q1263 . a. C 

Prandtl No. 

6,746 - 6.4478 

6. 8 - 6. 579 

6.55 - 6.308 

6. 8 - 6, 53 

Modified Grashaf. No+ = f t = 1* ), 

L 431 + 16 

1.37 + 15 

1. 26 + 14 

8.423 + 14 

Modified Gra&hof No. = ~ 

9.334 + 15 

8.978 + 14 

7.431 + 13 

6. 007 + 14 

Interface No. , I 

0.00756 

0.0218 

0.03118 

0.03633 : 

ULLAGE: density (Ib/It^) 

0,451 - 0*594 

0.459 - 0.546 

0.451 - 0.516 

0.455 . 0. 55' 

mass (lb) 

0.015-0.019 

0. 12 - 0. 143 

0.473 - 0.541 

0. 119 * Oi Ul 

specilic heat (Btu/Ib • F) 

0. 1S8 - 0.16 

0.159-0.16 

0.1 SB - 0. 159 

0. 158 - 0. 16. 
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Table A- 1 

- ULLAGE PROPERTIES (Page 1 of 2) 


, ac tfi 

6 in. 8G #14 

6 In. BG iilOS 

6 in. BG #3 

6 In. SG #4 

6 in. 8G #13 

6 In. 27C fflS 

03 

1. 140 + 03 

a. 867 + 03 

1.267 +04 

7. 59B + 03 

1. 14 + 03 

2. 838 + 04 

92 

4. 058 + 02 

2.948 + 02 

3. 609 + 02 

4.058 +02 

4,058 + 02 

4.028 + 02 

- 93. 329 

94.335 - 93. 562 

94.394 - 93,77 5 

94. 3.64 - 93, 618 

94.364 - 93.491 

93.449 - 93. 66 . . 

94. 174 - 93.241 

- 21.425 

21.61 - 21. 433 

12. 356 ^ 12. 275 

21.616 - 21,445 

21.616 - 21.416 

21.407 - 21. 455 

12.327 - 12.205 

Q.2Z7 

0.225 - 0. 227 

0.225 - 0.227 

0. 225 - 0. 227 

0. 225 - 0, 228 

0. 228 - 0. 227 

0.226 - 0-228 

0,039 

0,04 - 0.039 : 

0,04 - 0.039 

0,04 - 0.039 

0.04 - 0.039 

0,039 

0.039 

- 1. 127 

1. 194 - 1, 13 

1. 199 - 1, 147 

1. 196 - U 134 

1.196 - 1. 124 

1. 12i - 1.137 

I. 18 - 1. 105 

96 - 0.012049 

0.01265 - 0.01207 

0.0127 - 0.012226 

0.012678 - 0. 012 U2 

0.012678 - 0.0 120 22 

0.01199 - 0,01214 

0.01253 - 0.01185 

6. 57S 

6. 805 - 6. 587 

6. 822 - 6. 644 

6, 813 - 6. 60145 

6.813 - 6.568 

6.557 - 6.6126 

6. 757 - 6. 504 

^ 14 

1.228 +14 

1.202 + 14 

1. 354 + 15 

8,262 + 14 

1.247+ 14 

1.394 + 15 

14 

6.999 + 13 

7, 193 + 13 

a, 385 + 14 

5, 319 + 14 

8,211 + 13 

9. 531 + 14 


0.242 

0.03106 

0.0218 

0.0363 

0.242 

0.00977 

- 0, 559 

0.455 - 0. 566 

0.464 - 0.741 

0.456 - 0,782 

0,458 . 0.902 

0,456 - 0.765 

0.461 - 0.536 

' O.OIB 

0,015 - 0.019 

0.061 - 0.097 

0,015-0.026 

0,015 - 0.03 

0.015 - 0.025 

0.06-0.07 

- 0. 16 

0. I5B - 0. 16 

0. 159 - 0. 162 

0.158-0.162 

0,159-0.163 

0. 15B - 0. 162 

0- 159 - 0, 16 

n. 1C fl32 

12 in. IG m 

12 in. IG #15 

. 12 in. IG #35 



18 In. IG #D 

f 03 

7.6 + 03 

I, 14 + 03 

8. 85 1 + 03 

1,268 +04 

7, 597 + 03 

2. 859 + 04 

t 02 

4. 053 + 02 

4,013+02 

Z. 835 + 02 

3. 606 +02 

4,053 +02 

4. 069 + 02 

- .92. 39 

94. 31 - 93, 34 

94. 293 - 93, 416 

94. 34 - 93. 7 

94. 29 - 93. 449 

94.306 - 93,299 

94. 326 - 93. 55 

- 96, 755 

172,83 - 171.06 

172.8 - 171. 19 

98,79 - 98. 122 

172.795 - m.254 

1.72. 825 - 170. 979 

333.377 - 330. 633 

- 0, 231 

0.225 - 0.228 

0.225 - 0. 228 

0.225 -0.227 

0.225 - 0.228 

0.225 - 0.228 

0. 225 - 0- 227 

- 0,038 

0.04 - 0.039 

0.039 

0.04 - 0. 039 

0.039 

0.039 

0,040 - 0.039 

-1,044 

1, 191 - 1.113 

1, 19 - 1. 118 

i. 194 - 1, 14 

1.19-1.121 

1.191-1.109 

1.193-1.129 

7 - 0.0113 

0,01263 - 0.01192 

0,01262 - 0.01197 

0. 01265 - 0,0I;U7 

0.01262 - 0.01199 

0. 01263 - Q. 01189 

0.01265 - 0.01206 

6. 308 

6.S - 6. 53 

6.792 - 6. 549 

6. a - 6, 623 

6,79 - 6. 557 

6.796 - 6. 519 

6. aoz - (i. 584 

14 

a. 422 + 14 

1.2S1 + 14 

1. 212 + 14 

1,386+ 15 

8. 473 + 14 

1.349 + 15 

+ 13 

6.007 + 14 

8,186+13 

7.733 + 13 

9,673+ 14 

6.305 + 13 

7. 335 + 14 

S 

0. 03633 

0, 2423 

0. 0311B 

D.ozia 

0.03633 

0.00966 

- 0. S16 : 

0.455 - 0, 559 

0.458 - 0.549 

0.453 - 0.827 

0.460 - 0.828 

0.458 - 1,037 

0.445 - 0.513 

- 0, 541 

0. 119 - 0. 146 

0,12- 0, 144 

0, 474 - 0. 866 

0.126-0.217 

0. 120 - 0. 272 

1. 572 - 1. 814 

-0.159 

0.158-0,16 

0,159 - 0. 16 

0.158-0.163 

0.159 - 0.163 

0, 159-0.164 

0. 158 - 0. 159 
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Table 

FREON 113 PGA LIQUID - ULLAld 


Test No. 

6 in. 8G HISS 

6 In. BG !; 12S 

6 in. 8G i/5 

6 in, EC #16S 

Characteristic Temporature fCHTMP) ("F) 

8. B8 + 03 

1,269 + 04 

7.592 + 03 

8.868+03 

Characteristic Pressure fCHPRESS) (Ibf/in^) 

4.212 +02 

3, 609 + 02 

3.968 + 02 

2.948 +02 

LIQUID: 

density (Ib/ft^l 

94. 314 - 93,48 

94.282 - 93,556 

94.411 - 93.556 

94.392 - 93.749 


mass (IbJ 

12, 346 - 12.23t 

21. 598 - 21. 431 

21.627 - 21.431 

12.356 - 12.272 


spcciHc heat (Btu/lb *r} 

0.225 - 0.228 

0.225 . 0. 227 

0.225 - 0.227 

0.225 - 0.227 


thermal conductivity (Btu/hr ft *F1 

0.04 - 0.039 

0.039 

0.04 - 0.039 

O.04 - 0.039 


viscosity (Ib/ft hr) 

1.192 - 1. 123 

1. 189 - 1. 129 

1.201 - 1. 129 

1. 199 - 1. 144 


kinematic viscosity (ft^/hr) 

0.01264 - 0.0L2D1 

0,012613 - 0,012068 

0.012716 - 0.D12D68 

0.0127 - 0.01221 


Prandtl No. 

6,798 - 6, 565 

6.79 - 6. 585 

6. 827 - 6. 585 

6a 82Z - 6. 6^6 


Modified Grashof No, = f (qj^, = !• ). 

1. 253 + 14 

1. 367 + 15 

6. IH5 + 14 

1,207+ 14 


Modified Grashof No, r=l.}, Gr| 

7,069 + 13 

8.235 + 14 

5. 279 + 14 

7, 5D6 + 13 1 


Interface No, , 1 

O.OU06 

0,0218 

0,0363 

0,03106 

ULLAGE: 

density 

0. 46 - 0. 547 

0. 458 - 0- 564 

0.451 - 0. 568 

D. 456 - U, 745 


mass (lb) 

0,06 - 0.072 

0.015 - 0,018 

0.015 - 0.019 

0.06 - 0.098 


specific heat (Btu/lh ° F) 

0. 139 - 0. 16 

0. 159 - 0, 16 

0. 158 - 0. 16 

0. 158 . 0. 162 

Tqflt No. 

12 in. IG ff3 3 

12 in. IG i/30 

12 in. IG i/20 

12 in. LG i/26 

Characteristic Temperature (CHTMP) (°F) 

8. 859 + 03 

1.267 + 04 

7, 599 + 03 

1. 14 +03 

Characteristic Pressure (CHPRESS] 

4. 094 + 02 

3. 606 + 02 

4*nS3 + oz 

4,011 + 02 

LIQUID: 

density (Ib/ft^l 

94,282 - 93.426 

94.309 - 93.517 

94. 29 - 93. 364 

94.295 - 93.382 


mass (lb) 

98,732 - 97. 835 

172.83 - 171.378 

172,79'- - 171,099 

172.805 - 171. 132 


specific heat (Btn/ib °F) 

0.225 - 0.228 

6. 225 - 0. 227 

0.225 - 0.228 

0.225 - 0.228 


* icrmal conductivity (Btu/hr ft ”F) 

0-039 

0.040 - 6.039 

0.039 

0.039 


viscosity (Ib/ft hr) 

1,189 - UI19 

1, 191 - 1. 126 

1. 19 - i. 114 

1. 19 - 1. 115 


. kinematic viscosity (ft^/hr) 

0.01261 - 0.01254 

0.01263 - 0.01204 

0.01262 - 0.01193 

p. 01262 0.01174 


Prandtl No. 

6. 789 - 6. 551 

6. 797 - 6. 575 

6.791 - 6, 535 

6. 793 - 6, 54 


Modified Grnshol No. = f (q^> t = D, Gr^ 

1.259 + 14 

1.373 + 15 

8.401 + 14 

I. 257 + 14 


Modified Grashof No. = f Iq^, t = L), Gr^ 

7. 511 + 13 

9. 024 + 14 

5.744 + 14 

8, 563 + 13 


Interface No. , 1 

0. 03118 

O.OZiB 

0.03633 

0. 242 3 

ULLAGE; 

density (Ib/ft^l 

0.46 - 0. 520 

0.46 - 0. 547 

0.454 - 0-.554 

0.454-0,545 


mass (Ib) 

0. 481 - 0. 544 

0. 121 - 0. 143 

0.119-0.145 

0, 119 - 0, 143 


speciilc heat (Btu/lb ’F). 

0. 159 

6. 159 - 0. 160 

0:158-0. 16 

0. 158 - 0. 16 
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Table A-1 

L LIQUID - ULLAGE PROPERTIES (Page 2 o£ 2) 



6 in, 8G ffl6S 

6 in. SG 8llS 

6 in. SC 96 

6 in. 8G 88S 

6 ijj, IG fll 

6 in. IG #2 



8.068 + 03 

1,2654 + 04 

7. 603 + 03 

7. 581 + 03 

7. 596 + 03 

7.605 +03 



2.948 + 02 

3.609 + 02 

4.112 + 02 

1.968 + 02 

4. 004 + 02 

4.004 ^02 



94+ J92 - 9% 749 

94. 457 - 93, 692 

94.327 - 93, 394 

94,493 - 93. 651 

94. 38 - 93, 544 

94. 315 - 93.491 


.4^1 

12, 196-12.272 

21.638 - 21.463 

21. 608 - 21, 394 

21,646 - 21.453 

21. 62 - 21.428 

21.6-21.416 


n 

0.225 - 0. 227 

0. 229 - 0. 227 

0.225 - 0.228 

0.225 - 0.227 

0.225 - 0,227 

• 

0.22^ - 0. 22S 


9 


0.04 - 0. 039 

0.04 - 0.039 

0,04 - 0.039 

0.04 - 0.039 

0.04 - 0.0^9 


29 

I. 199 - 1. 144 

1.209 - 1. 14 

1. 193 - 1. 116 

1.208 - 1. 117 

1. 198 - I, 128 

1. 192 * 1, 124 


0.012068 

0,0127 - 0^41221 

0+01Z75 - O.OIZ166 

0.01265 - 0.01195 

0.01278 - 0.01214 

0.01269 - 0.01206 

0.01264- 0.01202 


89 

&4 822 - 6. 6^6 

6. 84 - 6. 621 

6. 802 - 6. 543 

6.852-6.61 

6.018 - 6. 582 

6, 799 - 6. 568 



1. 207 + 14 

1. 339 + 15 

8. 373 + 14 

8,073 + 14 

1.025 + 14 

1,033 + 14 



7. 506 + 13 

8.493 + 14 

5. 687 + 14 

5, 13 + 14 

6.415 + 13 

6. 36B + 13 



0.01)06 

0.0218 

0.0363 

0.0363 

0. 0363 

0.0363 


68 

0. 456 - 0, 745 

0.456 - 0.767 

0,450 - 0.964 

0.453 - 0. 554 

0. 449 - 0. 54 

0,455 - 0. 56l 


19 

0,06 - 0. 090 

0.015 - 0.025 

0.015 - 0.032 

0,015 - O.OlS 

0.315 - 0.018 

0. 015 - 0. 018 


6 

0,198-0.162 

0. 158 - 0. 162 

0. 159 - 0. 164 

0. 158 - 0. 16 

0. 158 - 0. 16 

0. 158 - 0. 16 


; #20 

I2in, lG(f26 

12 in, 1C #34 

12 in. IG 031 

12 in. IG 923 

12 in. IG #27 

18 in. IG 

18 in. IG #IS 


t. 14 + 03. 

8.858 +03 . 

1.267 + 04 

7, 594 + 03 

1 14 + 03 

2. 863 + 04 

3. 668 + 04 


4.013 + 02 

2. 8^1? 4 02 

3. 606 + 02 

4.053 + 02 

4+ 0 L3 ^ 02 

4,069 + 02 

5.214 + 02 

364 

94,295 - 93.382 

94.287 - 93, 613 

93.462 - 92.58 

94. 313 - 93.215 

94.276 - 93.236 

94. 253 - 93. 493 

94.261 - 93.218 

71. 099 

172.805 - 171. 132 

98,738 - 98.032 

171.278 - 169. 662 

172. 874 - 170-826 

172.769 - 170- 863 

333.117,330.432 

583. 006 - 376. 556 

28 

0.225 ■ 0,228 

0. 225 - 0. 227 

0, 228 - 0. 23 

0. 225 - 0. 228 

0. 225 - 0. 228 

0.225 - 0.228 

'.225 - 0.228 


0,039 

0, 039 

0.039 - 0.038 

0.04 - 0.039 

0,039 

0.039 

0.039 

4 

1.19-1,115 

1. 19 - 1, 134 

1. 122 - 1,057 

1. J94 - 1. 103 

1, 189 - 1. 104 

1. 187 - 1, 124 

1. 187 - 1. 103 

.01193 

0.01262 - 0.01194 

0,01262 - 0. 012109 

0.012 - 0.01142 

0.01265 - 0.01183 

0.01261 - 0.01184 

0.01259 - 0. 01202 

0.0126-0.01183 

39 

6.793 - 6, 54 

6, 79 - 6,6 

6. 561 - 6. 349 

6.804 - 6.498 

6.787 - 6. 503 

0.78 - 6,569 

6. 782 - 6, 498 


1. 257 + 14 

1, 227 + 14 

1. 383 + 15 

8. 565 + 14 

1,282 + 14 

1. 36 + 15 

1.395+ 16 


8.563 + 13 

B. 319 + 13 

9. 6B6 +14 

7. 075 + 13 

9:954 + 13 

7.238 + 14 

8, 357 + 15 


0. 2423 

0.03118 

0,0^18 

0. 03633 

0.242 3 

0,00966 

0,00754 

54 

0.454. 0,545 

0. 463 - 0. 844 

0.462 - Q. 9 

0,458 - 1,007 

0.463 - 0,812 

0.46 - 0, 531 

0. 465 - 0. 607 

45 

0, 119 - 0. 143 

0.4BS-O..8B4 

0, 121 - 0. 236 

0.12-0.264 

0.121-0.213 

1,624 - 1-.27S 

0.411-0.536 

6 

0,158-0.16 

0. 159 - 0.163 

0. 159 - 0. 163 

0. 159 - 0. 164 

0. 15° - 0. 163 

0,159-0.16 

0,159-0.16 























Appendix B 

TESt CONDITIONS AND TABULATED SENSIBLE HEAT DATA 

Tliis appendix contains th.e input heater heat flux, liquid {Freon. 113 PGA) 
level, and stratification test time for each test investigated. Included are 
the tabulated sensible heat gains by the test liquid, ullage, and the 
test tank structure^, the tank -wall, the foam insulation, and the silicone 
rubber heaters. This information is shown in Table B-1. 
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6 in. 8G 

6 in. 8G 

6 in. 8G 


Test No, 

No. 9S 

No. 2 

No. 1 


Liquid Level (percent) 

50. 0 

87. 5 

87.5 

1 

Heating Condition 

Liquid 

Liquid 

Liquid 


2 

Liquid Heater Electrical Source ql\ {Btu/hr ft ) 
Ullage Heater Electrical Source q^'^(Btu/hr ft^) 
Stratification Test Time (minute) 

701. 569 

tool. 99 

600. 85 

1 

1 

0. 0 

0. 0 

0. 0 

1 

3. 33 

2. 0 

3, 666 


Liquid Total Energy Supplied (Btu) 

45. 913 

65. 58 

71.966 


Ullage Total Energy Supplied (Btu) 

0. 0 

0. 0 

0.0 

' 1 

Summation — Total Energy Supplied (Btu) 

45. 913 

65. 58 

71.966 

■) 

Liquid — Measured Heat Flux qV (Btu/hr ft^) 

368. Z48 

582. 284 

389. 008 

J 

C 

Ullage — Measured Heat Flux q^* (Btu/hr ft*) 

1. 744 

2. 566 

1. 604 

■ . ,1 

Tank — Stainless Wall — Measured Sensible Heat (Btu) 

6. 089 

7. 867 

7. 244 

■ i 

Silicone Heaters — Measured Sensible Heat (Btu) 

4. 384 

5. 89 

5. 42 

1 

Foam Insulation — Measured Sensible Heat (Btu) 

8.675 

9. 095 

9. 233 


Freon 113 (PCA — Measured Sensible Heat (Btu) 

24. 102 

38, 1 1 

46. 669 

1 

Ullage — Measured Sensible Heat (Btu) 

0. 114 

0. 034 

0, 038 


Summation — Measured Absorbed Sensible Heat (Btu) 

43. 364 

60. 996 

68.604 

( 

1 

Percent — Measured Absorbed Heat/Total Energy Supplied 

94. 44 

93. 01 

95. 32 

1 


12 in. IG 

tZ in. IG 

12 in. IG 

12 

Test No. 

No. 25 

No, 32 

No. 22 

N.; 

Liquid Level (percent) 

87. 5 

50. 0 

87.5 

£ 

Heating Condition 

Liquid 

Liquid 

Liquid 

L* 

' 

Liquid Heater Electrical Source qU (Btu/hr ft^) 
Ullage Heater Electrical Source q”‘^(Btu/hr ft2) 
Stratification Test Time (minute) 

500. 56 

349. 77 

300. 077 


0. 0 

0, 0 

0. 0 

: 

. 8. 0 

13. 0 

15.0 

c 

Liquid Total Energy Supplied (Btu) 

524, 188 

357. 123 

589. 2 

5£ 

Ullage Total Energy Supplied (Btu) 

0. 0 

0. 0 

0. 0 


Summation — Total Energv* Supplied (Btu) 

524, 188 

357. 123 

589. 2 

5* 

Liquid — Measured Heat Flux q" (Btu/hr ft^) 
Ullage — Measured Heat Flux q”^ (Btu/hr ft^) 

328. 14 

206. 34 

214. 021 

2 

1. 194 

0. 723 

0.767 


Tank — Stainless Wall — Measured Sensible Heat (Btu) 

64, 858 

47. 752 

61.65 


Silicone Heaters — Measured Sensible Heat (Btu) 

23..58 

16. 912 

22.237 

: 

Foam Insulation — Measured Sensible Heat (Btu) 

72. 78 

67. 448 

72. 869 

. t 

Freon 113 PCA — Measured Sensible Heat (Btu) 

343.629 

210. 677 

420. 229 

3£ 

Ullage — Measured Sensible Heat (Btu) 

0. 25 

0.738 

0. 301 


Summation — Measured Absorbed Sensible Heat (Btu) 

505. 097 

342. 789 

577. 286 

5( 

Percent — Measured Absorbed Heat/Total Energy Supplied 

96, 35 

95. 98 

97. 97 

£ 
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Table B-1 

TEST CONDITIONS - SUPPLid 
SENSIBLE HEAT GAINS 


Table B- 1 

TIONS - SUPPLIED HEAT FLUXES - 
3LE HEAT GAINS (Page 1 of 2) 



6 in. 8G 
No. 1 

6 in. 8G 
No. 14 

6 in. 8G 
No. lOS 

6 in. 8G 
No. 3 

6 in. 8G 
No. 4 

6 in. 8G 
No. 13 

6 in. 27C 
No, IS 

6 in. 27G 
No. 3S 


87.5 

87.5 

50. 0 

87, 5 

97.5 

87.5 

50. 0 

87.5 


Liquid 

Liquid 

Liquid and 

Liquid and 

Liquid and 

Liquid and 

Liquid 

Liquid 




Ullage 

Ullage 

Ullage 

Ullage 




600. 85 

90. 127 

701. 569 

1 oor. 99 

600. 85 

90. 127 

2236. 686 

2884. 422 


0. 0 

0. 0 

701. 569 

1001. 99 

600. 85 

90. 127 

0. 0 

0.0 


3. 666 

25. 0 

2. 333 

2. 0 

3. 75 

25. 0 

1. 0 

1. 0 


71, 966 

73. 735 

32. 142 

65. 58 

73, 735 

73.735 

43.917 

94.392 


0.0 

0. 0 

32. 142 

13. 116 

1 4, 747 

14. 747 

0, 0 

0. 0 


71. 966 

73. 735 

64. 284 

78. 696 

88. 482 

88. 482 

43. 917 

94. 392 


389. 008 

51. 382 

419. 734 

620. 367 

386. 8l 

59. 366 

1529, 462 

1880. 114 

ii 

1,604 

0.249 

14. 206 

1 0. 772 

3. 338 

0, 809 

5.611 

8. 264 

7 ■ 

7. 244 

5. 421 

18. 8 

15. 867 

15. 06? 

8.756 

5. 91 

10. 355 


5.42 

3.892 

13. 396 

9. 65 

9. 049 

5.612 

4. 146 

7. 8 


9. 233 

9,116 

12. 54 

11.16 

1 1, 099 

10. 219 

8. 949 

10. 438 


46. 669 

42, 037 

19. 23 

40, 603 

47, 469 

48.569 

29. 031 

61.527 

t 

0. 038 

0. 041 

0. 465 

0. 141 

0. 206 

0. 1 32 

0. 1 07 

0. 054 

> 

68. 604 

60. 507 

64. 43 1 

77. 421 

89. 89 

73. 288 

48. 143 

90. 174 


95. 32 

82. 06 

100, 22 

98. 37 

93. 68 

82. 82 

1 09, 62 

95.53 

- 

12 in. IG 

12 in. IG 

12 in. IG 

12 i.'. IG 

1 2 in. IG 

18 in. IG 




No, 22 

No, 15 

No.. 35 

No. 17 

No. 28 . 

No. D 




87. 5 

87.5 

50. 0 

87.5 

87.5 

5 0. 0 




Liquid 

Liquid 

Liquid and 

Liquid and 

Liquid and 

Liquid 

. 





Ullage 

Ullage 

U! lage 





300, 077 

45.012 

349. 77 

500, 56 

300. 077 

753. 15 




0. 0 

j. 0 

350. 35 

500. 56 

300. 124 

0. 0 




15. 0 

99. 0 

9. 0 

8.0 

15.0 

9. 0 




589. 2 

583. h 

247. 239 

524. 188 

589, 2 

1197. 835 




0. 0 

0.0 

247.648 

1 04. 84 

117.976 

0, 0 



1 

r 

589. 2 

583. 31 

494. 887 

629. 028 

707. 176 

1197. 835 




214. 021 

29. 444 

223.179 

349. 47 

223. 311 

409.627 




0,767 

0. 101 

7, 55 

5. 763 

5. 844 

1. 639 



1. 

61, 65 

39. 912 

168.205 

. 109. 404 

125. 084 

201.998 



' 

22, 237 

14, 22 

60. 041 

34. 077 

37. 234 

49. 189 




72,869 

68.84 

105.571 

86.48 

89.537 

238. 016 




420. 229 

381.562 

157. 756 

365. 964 

438.47 

651.483 




0.301 

0. 262 

5.336 

1.207 

2. 295 

2, 607 




577.286 

504.796 

496. 909 

597. 132 

692, 62 

1 143. 293 




97.97 86.53 


94. 9Z 97.94 95.44 


100. 4 
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Talkie B-* 

TEST CONDITIONS - SUPPt| 
SENSIBLE HEAT GAINl 



6 in. 8G 

6 in. 8G 

6 in. BG 


Test No. 

No. 15S 

No. 12S 

No. 5 


Liquid Level (percent) 

50. 0 

87. 5 

DO 


Heating Condition 

Liquid 

Liquid 

Liquid 

. i 

2 

Liquid Heater Electrical Source qL (Btu/hr ft ) 
Ullage Heater Electrical Source q" (Btu/hr ft^) 
Stratification Test Time (minute) 

701. 57 

1001. 99 

600. 85 


0. 0 

0. 0 

0. 0 


3, 333 

2. 0 

3. 667 


Liquid Total Energy Supplied (Btu) 

45. 92 

65.58 

72.095 


Ullage Total Energy Supplied (Btu) 

0. 0 

0. 0 

0. 0 


Summation — Total Energy Supplied (Btu) 

45. 917 

65. 58 

72. 095 


Liquid — Measured Heat Flux q'' (Btu/hr ft^) 

395. 894 

603. 703 

387. 535 


Ullage — Measured Heat Flux qj^^ (Btu/hr ft^) 

1. 903 

2. 952 

1.8 


Tank — Stainless Wall — Measured Sensible Heat (Btu) 

6. 222 

8. 044 

7. 4£7 


Silicone Heaters — Measured Sensible Heat (Btu) 

4. 464 

5. 977 

5. 523 


Foam Insulation — Measured Sensible Heat (Btu) 

8. 827 

9. 267 

9. 263 


Freon 113 PCA — Measured Sensible Heat (Btu) 

25. 911 

39. 512 

46.5 


Ullage — Measured Sensible Heat (Btu) 

0. 125 

0. 039 

0.043 

" 3 

Summation — Measured Absorbed Sensible Heat (Btu) 

45. 549 

62. 839 

68. 796 

■ 1 

Percent — Measured Absorbed Heat/Total Energy Supplied 

99. 19 

95. 82 

95. 42 

1 


12 in. 1C 

12 in. IG 

1 2 in. 1 G 


Test No. 

No. 33 

No, 30 

No, 20 


Liquid Level (percent) 

50., 0 

87. 5 

87. 5 


Heating Condition 

Liquid 

Liquid 

Liquid 


2 

Liquid Heater Electrical Source (Btu/hr ft ) 

Ullage Heater Electrical Source q" (Btu/hr ft^) 
Stratification Test Time (minute) 

349. 77 

500. 56 

300, 077 


0. 0 

0. 0 

0. 0 


13, 0 

8. 0 

15. 0 


Liquid Total Energy Supplied (Btu) 

357. 123 

524. 188 

589. 2 


Ullage Total Energy Supplied (Btu) 

0. 0 

0. 0 

0. 0 


Summation — Total Energy Supplied (Btu) 

357. 123 

524. 188 

589. 2 


Liquid — Measured Heat Flux qV (Btu/hr ft^) 

208. 707 

362. 662 

205. 171 


Ullage — Measured Heat Flux (Btu/hr ft^) 

0. 657 

1. 192 

0. 74 


Tank — Stainless Wall — Measured Sensible Heat (Btu) 

41. 338 

64. 858 

59. 513' . 


Silicone Heaters — Measured Sensible Heat (Btu) 

14. 729 

23. 568 

21, 484 


Foam Insulation — Measured Sensible Meat (Btu) 

66. 745 

73. 774 

77. 077 


Freon 1 13 PCA — Measured Sensible Heat (Btu) 

213, 093 

314. 63 

402. 851 


Ullage — Measured Sensible Heat (Btu) 

0.671 

0. 250 

0. 291 


Summation — Measured Absorbed Sensible Heat (Btu) 

336. 576 

477. 08 

561. 216 


Percent — Measured Absorbed Heat/Total Energy Supplied 

94. 24 

91. 01 

95.25 







Table E- 1 

TTIONS - supplied HEAT FLUXES - 
^LE HEAT GAINS (Page a of 2) 


3 

6 in, 8G 

6 in. 8G 

6 in. 8G 

6 in. 8G 

o in. 8G 

6 in. IG 

6 in. 1 G 


5 

No, .5 

No. 16S 

No. US 

No. 6 

No. 8S 

No. 1 

No. 2 



87.5 

50. 0 

87. 5 

87. 5 

87. 5 

87. 5 

87, 5 


1 

Liquid 

Liquid and 

Liquid and 

Liquid and 

Liquid and 

Liquid and 

Liquid and 




Ullage 

Ullage 

Ullage 

*OHa;teF 

^4Laig'««i 



9 

600. 85 

701.569 

1001. 99 

600. 85 

600. 35 

600. 85 

600. 85 



0. 0 

701.569 

1001. 99 

600. 85 

0. 0 

0. 0 

0. 0 



3.667 

2. 333 

2. 0 

3. 8 

3. 667 

3, 7 

3. 7 


3 

72. 096 

32. 142 

65 58 

74. 718 

72. 095 

72.752 

72. 752 



0.0 

32. 142 

13.116 

14,944 

0. 0 

0, 0 

0. 0 


3 

72. 095 

64. 284 

78. 696 

89. 662 

72. 095 

72. 752 

72. 752 


03 

387. 535 

436. 428 

635. 42 

408. 07 

381. 824 

375. 963 

370, 305 


32 

1. 8 

10. 66? 

10. 178 

9. 723 

1.517 

1. 356 

1.594 


14 

7.467 

19. 823 

] 5. 645 

15. 912 

6. 978 

8. 4 

8. 488 


77 

5. 523 

14. 182 

9. 535 

9. 607 

5. 167 

6. 266 

6. 233 


37 

9, 263 

12.896 

1 1. 26 

11. 375 

9. 031 

9. 633 

9. 758 


12 

46.5 

19. 995 

41. 588 

50. 745 

45. 815 

45. 522 

44. 837 


39 

0. 043 

0. 488 

0. 133 

0. 242 

0. 036 

0. 033 

0. 039 


39 

68. 796 

67. 384 

78, 161 

87, 881 

67.027 

69. 854 

69. 355 


1 

95.42 

104. 82 

99. 32 

98. 01 

92, 07 

96. 01 

95. 33 


IG 

12in. IG 

i 2 in. IG 

12 in. IG 

12 in. IG 

12 in. IG 

12 in. IG 

18 in. IG 

18 in. IG 

0 

Nn. 20 

No. 26 

No. 34 

No, 31 

No. 23 

No. 27 

No. BS 

No. IS 


87. 5 

87.5 

50. 0 

87. 5 

87. 5 ' 

87. 5 

50. 0 

87. 5 

L 

Liquid 

Liquid 

Liquid and 

Liquid and 

Liquid and 

Liquid and 

Liquid 

Liquid 




Ullage 

Ullage 

Ullage 

Ullage 



i) 

300. 077 

45. 012 

349. 772 

500. 563 

300. 077 

45. 0116 

753. 151 

965. 0435 


0, 0 

0. 0 

350. 33 

500. 563 

300.424 

45. 0637 

0. 0 

0. 0 


15. 0 

99. 0 

9. 0 

8. 0 

15. 0 

99. 0 

9. 0 

9, 0 

38 

589. 2 

583. 308 

247. 239 

524. 188 

589.2 

583. 308 

1197, 835 

2558. 06 


0. 0 

b. 0 

247. 239 

104. 838 

117.976 

116, 797 

0, 0 

0. 0 

38 

589. 2 

583. 308 

494. 887 

j29, 026 

707. 176 

700. 105 

1197. 835 

2558. 06 

32 

2'^5. 171 

30. 664 

237. 126 

350. 637 

247. 854 

34. 946 

400. 823 

578. 124 

92 

0. 74 

0.102 

7. 69 

7. 732 

5,459 

0. 472 

1. 723 

2. 73 

58 

59.513 

40. 625 

162. 86 

128, 648 

126.866 

73.411 

206. 696 

302, 998 

58 

21.484 

14. 514 

57, 597 

38. 45 

37.778 

22. 89 

49. 245 

74. 144 

74 

77. 077 

75. 495 

103.787 

90. 176 

90. 516 

77. 084 

243, 66 

276. 177 

3 

402. 85 1 

397, 382 

167.614 

367. 186 

486. 661 

45Z. 867 

637.481 

1532. 445 

50 

0.291 

0. 264 

5. 436 

1,619 

2. 144 

I. 223 

2. 740 

1, 447 

3 

561. 216 

528. 28 

497. 294 

626. 079 

743.965 

627.475 

1139. 822 

2187. 21 1 

1 

95.25 

■90.56 

100,48 

99. 53 

105. 2 

89. 62 

95.15 

85. 5 
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